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1 INTRODUCTION 



ABSTRACT 

Analysis of the globular cluster (GC) system of the giant elliptical (E3) galaxy NGC 
4365, from eight Hubble Space Telescope/ Advanced Camera for Surveys (ACS) point- 
ings and a wide-field Subaru/SuprimeCam (S-Cam) image, is presented. Using magni- 
tude, colour and size criteria we obtain a catalogue of GC candidates. We also measure 
the photometric properties of the galaxy starlight, including a new measure of the ef- 
fective radius. We determine the lower limit on the number of GCs to be 6450 ±110 
and show that the GC system extends beyond 134 kpc (9.5 galaxy effective radii). We 
revisit the question of whether NGC 4365 has a bimodal or trimodal GC colour distri- 
bution and find support for three distinct GC colour subpopulations (i.e. blue, green 
and red). Sersic profile fits to the radial surface density of each subpopulation reveal 
that the blue GCs are more extended than either red or green GCs. The median half 
light radii for GCs in the blue, green and red subpopulations are 4.1±g-2, 3.0±q-J and 
2.8±o J parsec respectively. The estimated subpopulation ellipticities are 0.66 ± 0.06, 
0.55 ±0.07 and 0.16 ±0.25 for the blue, green and red GCs, where alignment with the 
photometric position angle of the galaxy (~ 42°) is assumed. A KS test on the mass 
functions show a > 98 per cent probability that all three subpopulations are distinct 
from one another. We also find radial gradients of GC size and colour (metallicity) 
and a blue tilt. The properties, including surface density profile, position angle, ellip- 
ticity and radial colour gradient, of the red GC subpopulation are very similar to the 
properties of NGC 4365's starlight. This result supports the hypothesis that red GCs 
are formed along with the bulk of the diffuse starlight in the galaxy. NGC 4365 has 
a kinematically distinct core and a significant misalignment between the photometric 
and kinematic major axes. We discuss the possibility that these kinematic features 
are related to the presence of the distinct third GC subpopulation. We briefly discuss 
implications for the formation of NGC 4365, finding that major merger, multi-phase 
collapse and accretion formation scenarios could all account for the existence of a third 
GC subpopulation. 



1.1 Globular cluster colour bimodality 

Bimodality in the colour distribution of GC systems in ellip- 
tical galaxies was first clearly demonstrated by |Zepf fc Ash-| 
|man| ( |T993 [). Since then, it has been shown to be a common 



often reveal distinct kinematic distributions ( Lee et al.|2010 
[Arnold et al.||201l[ [Foster et al.||2011[) similar t o those ob- 
served in the Milky Way's GC system ( Zinn 1985 1. The mean 



feature in large galaxies of all Hubble types (e.g. Forbes, 



Brodie fc Grillmair| p^97l [Larsen et a"L][200T| [Peng et al. 
2006). Spectroscopic observations of small samples of extra- 



colour/metallicity of the two subpopulations correlate with 
the host galaxy luminosity ( jPeng et al. 112006] ). The pres- 
ence of distinct subpopulations indicates that GC system 
formation was not a simple one-step process with (at least) 
two phases of GC formation suggested. Various models have 



galactic GCs have generally found them to be very old (~ 
12 Gyrs), similar to their Milky Way counterparts. Such 
old ages support the interpretation of colours as a good 
proxy for metallicity. Thus the photometric observations of 
large samples of GCs indicate two metallicity subpopula- 
tions: a blue, or metal poor, subpopulation and a red, or 
metal rich, subpopulation. These metallicity subpopulations 



been proposed to explain the two formation modes ( Ashman 
fc Zep f T992 , Forbes ct al. 1997| |C6te et al.||1998| [Beasley 
et al.|2002 , Bckki ct al.,2008j ). 



The interpretation of colour bimodality correspond- 
ing to metallicity bimodality has been challenged recently 



by Yoon et al. (2006), Cantiello & Blakeslee (2007) and 



Blakcslee et al. (2010) who suggested that the transfor- 



mation from colour into metallicity is strongly non-linear. 
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Beyond the Milky Way, there are very few GC systems ob- 
served with reasonable numbers of GCs in which to examine 
the spectroscopic metallicity distribution directly. However 
for those that do exist, i.e. NGC 4472 (|Strader et al.|2007l), 



NGC 5128 (Woodley et al. 20101 and NGC 4594 (Alves- 



Brito et al. 20111 spectroscopic metallicity bimodality is 



confirmed. 

With large photometric samples one can also examine 
other properties of extragalactic GC subpopulations. These 
include their radial and azimuthal distributions, luminos- 
ity or mass functions and physical sizes (if imaging with 
suitable resolution is available). Such properties also often 
reveal differences between the blue and red subpopulations 
(see review by Brodie fc Strader|2006 1. For example, the red 
GCs are usually more centrally concentrated than the blue 
GCs, and have a density profile, mean colour and spatial 
distribution that matches the bulge/spheroid component of 
the galaxy ( [Forbes fc Forte|2001| . The blue GCs have been 
measured to have larger mean sizes on average than the red 
GCs ( |Jordan|[2004l ). 

Although most large galaxies with well-studied GC sys- 
tems reveal clear colour bimodality, there are some excep- 
tions (see selected galaxies in Peng et al.|2006 1. One notable 
exception is the giant elliptical galaxy NGC 4365 which re- 
veals a broad colour distribution at all GC magnitudes. Here 
we revisit the issue of a possible third GC subpopulation at 
intermediate colours in NGC 4365. 



1.2 GC subpopulations of NGC 4365 

Here we study NGC 4365, a giant elliptical (E3) galaxy ( |Fer- 



rarese et al. 



Mpc away 



2006 1 on the far edge of the Virgo Cluster, 23.1 
Blakeslee et al.||2009|) with Mb = -21.3 mag 



i jFerrarese et al. 2006). There is debate as to whether the 
apparent Kinematically Distinct Core (KDC) of the galaxy 
is an actual kinematic property or the observed effect of 
triaxial orbits (e.g. [Davies et al.|[200T] [van den Bosch et al.| 
2008 1 . NGC 4365 is also one of few galaxies that rotate about 



the photometric minor rather than major axis (Davies et al 
2001 ). Tal et al. ( 2009 1 claim evidence of a faint fan feature 



in the SW of the galaxy. 



Forbes ( 1996 1 and Forbes et al. ( 1996 1 studied the GC 



system properties of NGC 4365 and other giant elliptical 
galaxies with KDCs. Bimodality was not suspected at the 
time in most GC systems and they found no significant 
differences between NGC 4365 and other giant elliptical 
galaxies with KDCs. Two independent groups working with 



HST/WFPC2 im aging, [Larsen et al.| ( |200l| ) and [Kundu fc 
|Whitmore| ( [200T| ), found that the GC system of NGC 4365 
was best fit by a unimodal distribution, whereas most of 
the elliptical galaxies in their samples were best fit by bi- 
modal distributions. Soon afterwards, Puzia et al. ( 2002 1 
analysed the GC system of NGC 4365, using the combina- 
tion of HST/WFPC2 optical and VLT/ISAAC K-band pho- 
tometry to break the age-metallicity degeneracy, which is a 
problem inherent to photometric analysis of GCs. In their 
analysis they found what appeared to be a slightly younger 
and very metal rich subpopulation in addition to two old 



indices (obtained from Keck/LRIS) of 14 GCs in NGC 4365 
it still looked likely that the galaxy hosted some old GCs 
and some very metal rich, young GCs that "conspire to pro- 
duce the single broad distribution observed in optical colors" 
( Larsen et a"L||2003[ ). 



However, when Brodie et al. ( 2005 I extended the sam- 



ple of NGC 4365 GCs with Lick indices by adding 19 new 
objects (also re-observing 3 objects from Larsen et al.|2003 l 
they found a uniform old age for all the GCs in the sample 
and suggested that the third GC subpopulation in the sys- 
tem, also found at intermediate optical colours, was due to 
a population with intermediate metallicity not young age. 
Hempel & Kissler-Patig (2004) published "better age deter- 



minations" for the intermediate GC subpopulation by in- 
cluding observations in another infrared filter and |Kundu| 
et al. ( 2005 I approached the problem with deeper HST ob- 



servations in the H-band. Larsen et al. ( 2005 1 confirmed 



the presence of intermediate optical colour GCs at small 
galactocentric radii (from an analysis of optical HST/ ACS 
photometry) and compared photometric age determinations 
with spectroscopic ones to assess the accuracy of photomet- 
rically determined ages. This work cast some doubt on the 
accuracy with which ages and metallicities can be measured 
using a combination of infrared and optical photometry. 
Later, Hempel et al. ( 2007 1 published a photometric compar- 



ison of elliptical galaxies in the centres of clusters and those 
in smaller groups, claiming (from optical and near-infrared 
photometry) that many group elliptical galaxies (including 
NGC 4365) host intermediate age GCs. Using g, z and K 



filter photometry, Chies-Santos et al. (20111 showed that 



subpopulations. When Larsen et al. (20031 published Lick 



NGC 4365's GCs all have similar ages (much like the other 
large ellipticals in their sample) but that the distribution 
in the g-z direction was significantly different to other large 
ellipticals. 

While there is still debate on the nature (in age, colour 
and metallicity) of a third GC subpopulation in NGC 4365, 
there is consensus that the GC system of NGC 4365 is dif- 
ferent to other galaxies of similar luminosity. Most analyses 
of NGC 4365's GC system subpopulations to date have been 
done with tens or hundreds of GCs. We revisit this issue in 
colour using thousands of GC candidates from optical pho- 
tometry over most of the spatial extent of the galaxy, to 
analyse whether there is statistically significant indication 
of three colour subpopulations in its GC system. 



1.3 Paper Structure 

In Section 2 and 3 we overview the data acquisition and re- 
duction for the Subaru/S-Cam and HST/ACS photometry 
respectively. We describe the GC candidate selection criteria 
and detail the photometric properties of the GC candidate 
samples in Section 4. Section 5 describes the properties of 
the galaxy light and we proceed to the analysis of the GC 
system properties in Section 6, first in terms of the properties 
of the total GC system and then in terms of the subpopula- 
tions and their characteristics. We discuss results in Section 
7 before concluding. 
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Figure 1. The central 34 X 25 arcmin section of the Subaru/S-Cam g' , r' and i' filter combined image (trimmed from 35 X 27 arcmin for 
cosmetic purposes) and the footprints of the 8 HST/ACS pointings are shown. The scale of the ACS pointings is shown in the bottom 
left corner. Image is centered on a=12:24:26.824; (5=+07;19;03.52 (J2000.0). At a distance of 23.1 ± 0.8 Mpc ( |BlakesIee et al.||2009| 
larcsec = 0.112 kpc. 



2 SUBARU/SUPRIME CAM DATA 

2.1 Observations and data reduction 

On 2008, April 1"* we obtained 35 x 27 arcmin tliree-filter 
imaging of NGC 4365 using Subaru/S-Cam. Conditions were 
good and the worst seeing was ~ 0.8 arcsec, exposure times 
were 5 x 130 s, 5 x 70 s and 5 x 60 s for g' , r' and i' filters 
respectively. The images were bias subtracted, flat field cor- 



rected and stacked using the SDFRED package ( Yagi et al, 
Ouchi et ar]|2004|| and put onto the USN0-B2 astro- 



2002 



metric system. Figure [ij shows a combination of the g' , r 
and i' filter images. The S-Cam instrument has a pixel scale 
of 0.202 arcsec. 

The NGC 4365 galaxy light was modelled and sub- 
tracted in each of the three filters using the IRAF task 
ELLIPSE before an object detection algorithm was em- 
ployed. This was done to increase the success with which the 
DAOFIND detection algorithm finds faint objects in the cen- 
tral regions of the galaxy (it is unsuccessful at finding objects 
in areas in which the background surface brightness varies). 
The standard deviation of the background, after galaxy sub- 
traction, was a = 8.64 counts in g' , a = 11.16 counts in r' 



and a = 13.15 counts in i' . Detection thresholds of 2.7cr, 
3.0(7 and 2.1a in g' , r' and i' respectively were used so as 
to probe as deep as possible in all areas of the images, with 
the assurance that the selection criteria employed later (see 
Section 4.2) would remove most spurious detections (due to 
cosmic rays, background fiuctuations and galaxy subtraction 
artifacts). 



2.2 Photometry and calibration 

The point- like objects (see Section 4.2.1 for details on point 
source determination) brighter than i' = 22 were used to 
define an optimum aperture for which to extract the pho- 
tometry. The photometry was extracted at a radius of 3.5 
pixels in g' and 3.0 pixels in r' and i' (the seeing was ~ 0.6 
in r' and i' but ~ 0.8 in the g' filter). The radius of extrac- 
tion was chosen as a trade-off between increasing error from 
background noise and increasing uncertainty in the correc- 
tion for light outside the extraction aperture (i.e. an aperture 
correction). The aperture correction was — 0.382±0.016 mag 
in g\ -0.263 ± 0.012 mag in r' and -0.249 ± 0.011 mag in 
i' . The sky value was calculated as the mode of all pixels in 
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Figure 2. Calibration of S-Cam pliotometry to SDSS filters. Tfie difference between zeropoint corrected S-Cam and SDSS magnitudes 
for g' , r' and i' filters are plotted from top to bottom. Left: The residuals are plotted against the {g' — i') colour and the colour range 
is restricted to that expected for GCs. All three filters show asymmetry in their distributions, scattering to brighter S-Cam magnitudes. 
Right: The residuals are shown against SDSS magnitude. 



an annulus between 15 and 20 pixels from the centre of the 
object. 

Standard stars were not observed with NGC 4365, in- 
stead the photometry was calibrated using bright point-like 
objects (18 < i' < 22) detected in the NGC 4365 field that 
were also found in the Sloan Digital Sky Survey (SDSS) cat- 
alogue. Both PSF and model derived SDSS magnitudes were 
used for this and no systematic difference was found in the 
results. There were 1822, 2174 and 2813 cross matched ob- 
jects for calibration in g' , r' and i' filters respectively and the 
photometric calibration zeropoints were determined from a 
best fit linear relation between the SDSS catalogue magni- 
tudes and the S-Cam instrumental magnitudes. These were 
found to be zpg, = 27.64 ± 0.03, zpr> = 27.76 ± 0.08 and 
zpi' = 27.72 ± 0.10 on the AB photometric system. The ze- 
ropoint corrected residuals of point-like objects in each filter 
are shown in Figure [2] plotted against g' — i' colour (left) 
and filter magnitude (right). The SDSS and S-Cam g' , r' and 
i' filters show good agreement with no obvious systematic 
colour or magnitude trend. 

It was necessary to correct for foreground dust extinc- 
tion as such a correction was not applied to the SDSS cal- 
ibration objects. We used the [Schlegel et al.| ( |1998| ) dust 
maps to calculate the extinction correction in each filter. 
The values were compared at different positions across the 
field and not found to vary significantly. For g' , r' and i' fil- 
ters the extinction correction values used were Agi = 0.081, 
Ar' — 0.060 and Aii — 0.045 mag. Hereafter we quote ex- 
tinction corrected magnitudes and colours. 



3 HST/ACS DATA 

3.1 Observations and object measurement 

We obtained, from the Hubble Legacy Archive, g and z {g 
and i for one pointing) filter imaging for eight separate point- 



Table 1. Imaging obtained from the Hubble Legacy Archive. 





Central 
R.A. Dec. 


Exp 
Filter 


osures 
Time (s) 


HST 
ID 


1 


12:24:27.0 


07:19:20.8 


F475W 
F850LP 


750 
1210 


9401 


2 


12:24:17.5 


07:21:09.5 


F475W 
F850LP 


680 
1170 


10582 


3 


12:24:16.0 


07:17:37.1 


F475W 
F850LP 


680 
1170 


10582 


4 


12:24:26.3 


07:15:44.8 


F475W 
F850LP 


680 
1170 


10582 


5 


12:24:36.4 


07:17:25.1 


F475W 
F850LP 


680 
1170 


10582 


6 


12:24:41.3 


07:20:18.2 


F475W 
F850LP 


680 
1170 


10582 


7 


12:24:31.5 


07:22:20.4 


F475W 
F850LP 


680 
1170 


10582 


8 


12:24:07.2 


07:12:11.3 


F475W 
F775W 


1744 
1624 


9488 



ings of the HST/ACS instrument around NGC 4365. Table[T] 
summarises these observations and Figure [l] shows the ACS 
pointing footprint on the S-Cam image. These archival data 
probe down to z = 25.2, at a 50 per cent completeness level 



( Jordan et al.|2007[ ). The resolution of HST/ACS data (pixel 
scale is 0.05 arcsec) partially resolves GCs at the distance 
of NGC 4365, 23.1 Mpc away ( [Blakeslee et al.|2009| ). 

The eight individual ACS fields were analysed using a 
custom built pipeline (see e.g. Strader et al.,2006, and . Spitler] 
et al.|[2006 l to find small, round objects and measure their 
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Figure 3. Size in parsec vs. g—z colour distribution for GC candi- 
dates brighter than the turnover magnitude of z = 23.4, detected 
in the HST/ACS fields. The GC candidate selection criteria is 
marked with a box. The blue GC candidates tend to have larger 
average sizes and the red GC candidates have a larger dispersion 
in colour. 



magnitudes and half light radii (henceforward referred to as 
object size). For details on the methods used by the pipeline 
including point spread function determination see |Strader| 



et al. (2006 1. The object lists from each field were collated, 



with the arithmetic mean taken of the sizes and magnitudes 
in the field overlap areas, and associated errors adjusted to 
refiect the more accurate measurement. 

The measured magnitudes and sizes are compared with 



those published for objects in the central field (Jordan et al 



20091. The magnitude measurements show no evidence of a 



statistically-significant systematic offset from the published 
data. The scatter increases to fainter objects but is on the 
order of 0.1 magnitudes. The size measurements in this work 
show a small ~ 0.005 arcsec systematic offset to smaller sizes 
attributed to differences in size measurement techniques. 



4 GC CANDIDATE SELECTION 
4.1 HST/ACS GCs 

GC candidates were selected from the HST/ACS images 
based on their size, colour and magnitude. The GC candi- 
date distribution is selected to have colours 0.7 < g — z < 1.6 
and sizes 0.1 pc < r?i < 20 pc. The choice of g — z colour up- 
per and lower bounds was based on the colour size diagram 
shown in Figure [3] There is a clear drop off in the density of 
objects bluewards of 0.75 and a similar but less clear drop 
off redwards of 1.6 (-2. 4 < [Fe/H] < 0.1 8 using the em- 
pirical transformation of Peng et al.| ( |2006| ). While most of 
the GC candidates in Figure [3| are resolved, the redder GC 
candidates show a size distribution that overlaps with ob- 
jects of zero size (objects indistinguishable from stellar point 
sources). A size cut at 0.1 pc (10~^ arcsec) was made because 
at smaller sizes the objects did not show a distribution that 
was centrally concentrated on the galaxy and therefore does 
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Figure 4. Colour magnitude diagram for GC candidates selected 
from HST/ACS imaging. 



not sample mainly GCs. A generous upper bound on size 
of 20 pc (0.2 arcsec) was employed to exclude background 
galaxies but include possible Ultra Compact Dwarfs (UCDs) 
in the system. Lastly, only objects fainter than z — 19 
(M2 = —12.8 mag) are considered GC candidates of NGC 
4365, following the convention of the ACS VCS published 



catalogue (Jordan et al. 20091. The catalogue contains no 



objects fainter than z = 26.2 mag. Objects brighter than 
z = 19 are most likely to be stars. 



4.2 Subaru/SuprimeCam GCs 

4-2.1 Point source determination 

At the distance of NGC 4365, GCs are not resolved by 
our Subaru observations and consequently they cannot be 
separated from stars in our Galaxy via a size distinction. 
By selecting only objects that are unresolved in the S-Cam 
imaging background galaxies are excluded from the analysis 
of NGC 4365's GC system. Here the distinction between a 
point-like object (mostly stars and NGC 4365 GCs) and an 
extended object (galaxies) is determined by a measure of the 
flux difference between two aperture radii. Shown in Figure 
[5] are the objects that were identified to be point-like in all 
three filters. Magnitudes were extracted for two apertures 
different in radius by a half or full pixel (filter dependent) 
both centered on the object. Extended objects have extra 
light in the larger aperture compared to point-like objects. 

The lists of point-like objects for each filter were cross- 
matched in position and only objects classed as point-like in 
all three filters were kept in the sample. This was done us- 
ing the IRAF task TMATCH with a tolerance in positional 
offset between images of 1.26 arcsec (this was determined by 
the uncertainty in astrometry at the image edges). Because 
the photometry is deeper in the i' filter than either r' or 
g' filters this procedure likely excluded genuine faint GCs 
from the analysis that were detected in i' but not in either 
of the other filters. Using a cross-matching technique spuri- 
ous detections (due to cosmic rays etc.) in any of the three 
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Figure 5. The difference in flux between two aperture radii (the 
smaller of which is the extraction aperture) vs. aperture corrected 
magnitude at the extraction aperture for all objects found by the 
DAOFIND detection algorithm. Objects determined to be point- 
like in all three S-Cam filters are shown in black dots, overlapping 
the total object detection in the g' filter (grey dots). The dot- 
ted line shows the cutofi' between point-like objects and extended 
objects in the g' filter. Bright point-like objects have aperture 
differences of ~ 0.1 magnitudes but fainter objects show signifi- 
cant scatter in measured magnitudes. A difference in magnitude 
between the inner and outer aperture of up to ~ 0.6 mag is re- 
quired to include fainter point-like objects. 



filters were eliminated and only objects witli piiotometry in 
all three filters selected. 



4-2.2 Colour-colour and magnitude selection 

As shown in Figure [s] the HST/ACS is able to partially re- 
solve GCs at the distance of NGC 4365 and can therefore 
directly distinguish between unresolved foreground stars, 
GCs and resolved background galaxies. We used the central 
HST/ACS pointing to determine where on a S-Cam r' — i' 
vs. g' — r' diagram the GC sequence lies by matching S-Cam 
point source detections with ACS GC candidates. A line was 
fitted to the brightest S-Cam objects found to be GC can- 
didates in the ACS field and used to make a box 2cr in each 
direction from the best fit GC sequence line. Objects were 
considered to be consistent with the GC colour-colour defi- 
nition if the standard error in their photometry placed them 
within the box (see Figure[6|. The selected distribution con- 
tains objects with colours 0.6 < g — i < 1.4, corresponding 
to metallicities of 
|& Hwang (2010 



-2.1 < [Fe/H] < 0.7 when the |Lee, Park 
empirical transformation is used. 
An upper limit on the magnitude of any GC candidate 
was set to i' — 19 {Mi/ — —12.8 mag) by the brightest 
GC found by HST/ACS in the central field of NGC 4365. 
A limit of 0.1 mag in photometric error resulted in a lower 
limit of i' m 25 mag on the blue end of the colour-colour 
selection and i' « 24.5 mag on the red end of the colour- 
colour selection. 
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Figure 6. Colour selection of GCs in the S-Cam field-of-view 
based on the r' — i' vs. g' — r' colour of matched objects in the 
HST/ACS data. The objects shown are all brighter than i' = 
22 mag and have magnitude errors smaller than 0.02 for visual 
clarity. S-Cam detected point-like objects are plotted in black 
plusses and HST/ACS matched objects are overplotted as red 
dots. The HST/ACS matched objects define a very tight sequence 
(black line) that is used to define the area from which S-Cam 
objects are selected as GC candidates (red box). 
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Figure 7. Colour magnitude diagram for GC candidates selected 
from S-Cam imaging. 



4.3 Contamination and completeness 

At galactocentric radii smaller than 0.9 arcmin S-Cam de- 
tection numbers drop due to galaxy subtraction artifacts 
and at radii greater than 3.4 arcmin ACS detections are 
no longer spatially complete due to the tiling pattern of 
HST/ACS pointings. Therefore only in an annulus between 
0.9 and 3.4 arcmin from the centre of NGC 4365 is it possi- 
ble to compare the object detections and photometry for the 
S-Cam and ACS imaging. This was done to determine the 
level of contamination in the S-Cam GC sample (contamina- 
tion is defined as non GC objects that meet the colour, size 
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Figure 8. Estimate of contamination in the S-Cam data by com- 
parison witli HST/ACS matched objects. The soUd line shows 
the number of S-Cam GCs in each magnitude bin, the dashed 
hne shows the number of those objects that are determined to be 
contaminants by HST/ ACS imaging and the dotted hne shows the 
percentage contamination of the S-Cam GC sample as a function 
of magnitude. 



Figure 9. Estimate of the completeness of the S-Cam photomet- 
ric sample by comparison with HST/ACS photometry. The solid 
line shows the number of HST/ACS GCs in each magnitude bin, 
the dashed line shows the number of those objects that were also 
detected in the S-Cam imaging and the dotted line shows the 
percentage completeness of the S-Cam photometry as a function 
of magnitude. 



and magnitude conditions set for the S-Cam photometry) 
as weU as to determine the photometric completeness of our 
S-Cam imaging. This analysis is relative to the HST/ACS 
sample but the ACS imaging is almost perfectly complete 
several magnitudes deeper than our Subaru/S-Cam imaging 
( Jordan et al.|2009| ). 

The contamination percentage of the sample in S-Cam 
was determined by keeping record of the GC candidates de- 
tected in the ACS imaging but removed by the conditions 
outlined in Section 4.1. At a magnitude of z = 24.2, 50 per 
cent of S-Cam GC candidates are determined to be contam- 
ination and at 2 = 24.5 mag essentially all of the S-Cam 
GC candidates are contaminants. Across the whole GC can- 
didate sample the contamination determined in this way is 
4.14 ± 0.35 objects arcmin^^ and if the sample is restricted 
to only objects brighter than z = 23.4 mag this value drops 
to 1.23 ± 0.19 objects arcmin"^. The luminosity function of 
GCs can be well described by a Gaussian distribution peaked 
at the turnover magnitude of ~ —8.4 mag ( [Jordan et al.| 



2007), which corresponds to z — 23.4 mag at the distance 



of NGC 4365 and all GC candidates described in further 
analysis are brighter than this turnover magnitude. 

The completeness of the S-Cam imaging was deter- 
mined by keeping record of how many of the HST/ACS 
detected objects were also detected by the Subaru/S-Cam 
imaging at different magnitudes. The S-Cam imaging is 
found to be 50 per cent complete at i' — 23.8 mag and 
70 per cent complete at i = 23.6 mag as shown in Figure |9] 

In the very central regions (radii < 0.325 arcmin) of 
NGC 4365 the HST/ACS detection of GC candidates is 
inhibited by the high surface brightness of the galaxy. In 



Jordan et al. (20071 the completeness of objects is tabulated 
as a function of object size, object magnitude and surface 
brightness of the galaxy at the location of the GC candi- 



date. For this analysis it is only necessary to consider com- 
pleteness as a function of radius. We measured the average 
surface brightness of the galaxy in 0.1 arcmin annuli on the 
ACS imaging. We compared that to the completeness frac- 
tion for an object of mean size = 0.039 arcsec (4.3 parsec) 
and z — 23.4 mag, which is representative of the GC can- 
didate sample. At distances of 0.2 and 0.325 arcmin from 
the galaxy centre the ACS imaging is 93.3 and 98.9 per cent 
complete. 



4.4 Overview of GC selection 

The selection of GC candidates from the Subaru/S-Cam 
photometry was done based on object size, magnitude and 
locus in colour-colour space compared to GC candidates ob- 
tained from HST/ACS photometry. Point-like objects were 
selected in each S-Cam filter using a limit on the flux dif- 
ference between two radii apertures (see Section 4.2.1). The 
astrometry of the ACS GC catalogue was shifted to that 
of the S-Cam catalogue. The positions of objects classed 
as point-like in all S-Cam filters were compared with the 
positions of GC candidates obtained from the ACS pho- 
tometry (see Section 4.1 for ACS selection criteria; colours 
0.7 < g — z < 1.6, sizes 0.1 pc < r^ < 20 pc and magnitudes 
26.2 < z < 19.0). The matched S-Cam point-like objects 
were used to determine the locus of GCs on a r' — i' vs. 
g' — r' diagram. No contaminating objects could be visually 
identified on inspection of the GC candidates brighter than 
i = 21 mag. Finally, estimates of the contamination and 
completeness of the S-Cam catalogue were determined by 
comparison with the ACS catalogue (see Section 4.3). The 
full catalogue of NGC 4365's GC candidates, found either 
ACS or S-Cam imaging, is included as online supplementary 
material in FullGlobList.dat. 
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Figure 10. Isophotal parameters for g',r' and i' filters from Subaru/S-cam photometry of NGC 4365. Prom left to right and top to 

bottom the parameters plotted are: surface brightness of the galaxy light in g' (erosses), r' (stars) and i' (plusses) filters; galaxy eolour; 
elliptieity of the fitted isophotes in all three filters as labelled before; position angle of the fitted isophotes in all three filters as labelled 
before; higher order Fourier coefficients to the sine (S3 and S4) and cosine (C3 and C4) terms. All parameters are plotted against 
equivalent galactocentric radius, calculated as the geometric mean of the semi-major and semi-minor axis lengths of the fitted ellipse. 
The fitted values become unstable beyond 3 arcmin. Isophotal parameters for the i' filter axe included as online supplementary material 
in GalLightTable.dat. 



© 0000 RAS, MNRAS 000, 000-000 



Wide-field imaging of NGC 4365's globular cluster system: The third subpopulation revisited 9 



Filter 




n 


Re 




(mag arcsec"-^) 




(arcmin) 


g 


21.87 ± 0.11 


6.02 ±0.21 


2.40 ±0.14 


r 


21.46 ±0.12 


5.92 ±0.23 


2.06 ±0.12 


i 


21.20± 0.11 


5.97 ± 0.24 


2.10 ±0.11 



Table 2. Values for g, r and i Scrsic profiles fitted to galaxy light. 

5 ANALYSIS OF THE GALAXY LIGHT 

5.1 Parameters of the isophotal fits 

We used the IRAF task ELLIPSE to model and subtract the 
galaxy hght from the Subaru/S-Cam photometry. From the 
model of the galaxy hght we obtain the surface brightness, 
ellipticity, position angle and higher order Fourier terms (S3, 
84, C3, C4) in g',r' and i' filters for NGC 4365. Results are 
plotted in Figure [lO] and briefly discussed here. We have 
surface brightness information for all three filters between 
0.1 and 5 arcmin. We fit Sersic profiles ( [Graham fc Driver] 



2005), i.e 



fJ,{R) = fle + 



2.5bn 

InlO 



(1) 



to the three surface brightness profiles, where J?e is the ef- 
fective radius of the galaxy, fie is the surface brightness at 
that radius, n is the shape parameter of the Sersic profile 
and 6„ ~ 1.9992n — 0.3271. The fitted values are tabulated 
in Tabled 

The quoted value for NGC 4365's effective radius 
(Goudfrooij et al. 1994 Bender et al. 19921 is 1.1 arcmin 



based on a de Vaucouleurs profile (equivalent to a Sersic pro- 



file with the n parameter set to 4) fitted by Burstein et al. 
(19871. When a de Vaucouleurs profile is fit to the data in 
this work the effective radius is found to be 1.32 ± 0.03 ar- 
cmin but we find this to be a significantly poorer fit than a 
general Sersic profile. For all further analysis we use our Re 
value of 2.1 ± 0.1 arcmin (14.1 ± 0.7 kpc), with an n value 
of 6.0 ± 0.2. This is close to the value of 1.6 ± 0.1 arcmin 
(10.9 ± 0.7 kpc) with n = 5 



(20101 



in z that Chen et al 
found and slightly smaller than the value of 3.07 ± 0.22 ar- 
cmin (20.6±1.5 kpc) with n — 7.1±0.4 in B that Kormendy 



et al. 



(2009| found. 

Results for the ellipticity, position angle and higher or- 
der Fourier terms in this work agree with the values |Goud-| 
al. (19941 found in the B,V and I filters and ex- 
tend 2 arcmin further, to ~ 1.57ie. The ellipticity measured 
in this work varies between 0.22 at 0.15 arcmin and 0.28 at 
2.5 arcmin and is consistent with that seen by [Goudfrooij | 
et al. (19941 interior to 1 arcmin. The position angle shows 
a small twist of < 5°beyond the radial range explored by 
Goudfrooij et al. (1994 1 . The S3 parameter (coefficient to 



the sine term in the fit) shows an indication of increase at 
the edge of their radial range (1 arcmin) which is confirmed 
in this work (between 1 and 2 arcmin) and correspondingly 
the S4 parameter increases between 1 and 2 arcmin. There 
is no significant difference in C3 and C4 parameters between 
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Figure 11. Surface density of GC candidates brighter than the 
turnover magnitude plotted against galactocentric radius. The S- 
Cam points (solid points) were calculated using objects brighter 
than i' = 23.6 and the ACS points (hollow diamonds) were cal- 
culated using objects brighter than z = 23.4. The innermost two 
ACS values were corrected for incompleteness using the values 
discussed in the text. No other normalisation correction has been 
made to either ACS or S-Cam radial surface density and the re- 
markable agreement in the region of overlap (between 1 and 3 
arcmin) is an indication that both samples are almost perfectly 
complete to the turnover magnitude and contain very little con- 
tamination. The line plotted is a fitted Sersic profile with a con- 
stant background value. See text for further details. 



parameter remains at values between -0.01 and -0.02 out to 2 
arcmin, indicating that NGC 4365 has boxy isophotal struc- 
ture to large radii. The galaxy colour becomes bluer with in- 
creasing galactocentric radius, changing by ~ 0.1 mag from 
~ 0.1 to 2.5 arcmin for both g ~ i and g ~ r colour indices. 
A similar trend, implying a negative metallicity gradient, is 



seen in the Goudfrooij et al. ( jl994| ) B-I colour index. 



Goudfrooij et al. (19941 and this work. The C4 disky/boxy 



6 ANALYSIS OF THE GC SYSTEM 
6.1 Characterising the total GC System 

6.1.1 Surface Density 

By combining the wide field imaging from S-Cam with the 
spatially resolved imaging from ACS and further restrict- 
ing the GC candidate sample to objects brighter than the 
turnover magnitude we can derive a very accurate radial 
surface density profile over a factor of 100 in radius. To de- 
rive surface density profiles, the number of GC candidates in 
each radial bin were divided by the total area in that radial 
bin. The results for both the S-Cam and ACS GC candidate 
samples are plotted on Figure [TTj 

We expect very little contamination in the density pro- 
file, with more than 90 per cent completeness in this mag- 
nitude range (see Section 4.3). A Sersic profile plus a back- 
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(arcmin~^) 



Re bg 

(arcmin) (arcmiri"^) 



3.9 ± 1.3 2.68 ±0.41 6.1 ± 1.2 0.89 ±0.13 



Table 3. Fitted values of the Scrsic fit to the GC surface density. 



ground parameter, i.e. 

P{R) = Pe exp ( -b„ 

where 



6„ = 1.9992n - 0.3271 



+ bg (2) 



(3) 



is fitted to the combined S-Cam and ACS radial surface 
density of GC candidates and shown in Figure Re is 
the effective radius of the GC system, is the density at 
the effective radius and n is the shape parameter of the 
Sersic profile. The fitted parameters are recorded in Table 
[3] The value for contamination expected from the analysis 
in Section 4.3 (1.23 ± 0.19 arcmin"^) agrees within 2a with 
the value from the Sersic profile fit (0.89 ± 0.13 arcmin^^). 

We are not able to define the edge of NGC 4365's GC 
system with certainty as the profile shown in Figure [Tl] is 
still decreasing at ~ 20 arcmin (134 kpc), which is the very 
edge of the spatial coverage available to us with S-Cam 
imaging. In order to confirm the edge of the GC system 
even larger field of view imaging would be required. We can 
make an estimate of the total number of GCs associated 
with NGC 4365 by simply doubling the number of GC can- 
didates brighter than the turnover magnitude (i' — 23.6 and 
z = 23.4). This assumes that there is no inherent asymmetry 
in the GC luminosity function, the GC system of NGC 4365 
contains few GCs beyond 20 arcmin and that the sample is 
complete and free of contamination. We determine the total 
number of GCs in NGC 4365's GC system to be 6450 ± 110. 
This is likely to be a lower limit and doubles the expected 
number of G Cs in NGC 4365 fr om the value of 3246 ± 598 
calculated by Peng et al. (20081. They use the narrow field 



of view of one HST/ACS and three HST/WFPC2 pointings 
to derive an extrapolated GC surface density profile that 
is not publicly available. We suggest that the surface den- 



sity profile that Peng et al. ( 2008 1 fit to the GC candidates 



underestimated the radial extent of NGC 4365's GC system. 

These total GC numbers convert to specific frequency 
I Harris fc van den Berghl|1981| ) values of Sn = 3.86 ± 0.71 
I Peng et al.|2008| ) and 3^= 7.75±0.13 determined from this 
work for My = —22.31. Several other large elliptical galaxies 
with similar My have specific frequency values comparable 
to the Sn measured here. NGC 14 07 has Sn = 7 98 ± 0.87, 
NGC 1399 has Sn = 6.72 ± 0.81 (ISpitler et al.pOOSl and 



Peng et al. (20081 find Sn = 5.2 ± 1.4 for M84. Fitting a 



power law to the radial density beyond the Sersic effective 
radius we find a slope (i.e. power-law exponent) of —1.21 ± 
0.03 to allow comparisons with analysis of other galaxies in 
the literature. 

By binning the GC candidate positions into R.A. and 
Dec. defined squares we derive a two-dimensional density 
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Figure 12. Surface density of GC candidates in two dimensions 
from S-Cam data. The colour ranges from red to blue indicating 
decreasing GC surface density. Ellipses showing the ellipticity and 
position angle of the GC system are plotted with dotted black 
lines and galaxy light isophotes are plotted for comparison in 
solid black lines. The GC system is much more elongated than 
the starlight. 



distribution for the GC candidates. This is shown in Figure 
T2I along with selected galaxy surface brightness isophotes 
from the i' filter as determined by the IRAF task ELLIPSE. 
McLaughlin et all (119941) use 



(j{R, e) = fc_R""[cos^(6' - PA) 

+ {l-e^y^sm^{e-PA)]-°'^^ +hg (4) 

to fit for the position angle {PA) and ellipticity (e) of the 
GC system of M87. They fit to the density of GCs as a 
function of azimuthal distribution and we follow the same 
procedure. We use the previously determined power law ex- 
ponent (q) and the fitted background value {bg), allowing 
the normalisation constant (k) as well as the position angle 
and ellipticity to vary. We find the position angle of the GC 
system to be 39.7±2.2° and the ellipticity to be 0.63±0.02, 
see Figure |13| The GC system of NGC 4365 has a similar 
position angle to the galaxy light (~ 42°) but is clearly much 
more elliptical than the starlight of NGC 4365 (which has 
an ellipticity of ~ 0.25). We see qualitative agreement with 
this analysis in Figure [T2I 



6.1.2 Colour Distribution 

Before performing a quantitative analysis of the surface den- 
sity characteristics of GC subpopulations in NGC 4365 we 
can make a qualitative assessment of colour-radius substruc- 
ture by examining the surface density of GCs in colour- 
radius space. This is shown in Figures [l4] and [l5] At galacto- 
centric radii larger than 2 arcmin in the top panel of Figure 
|14| the blue and red GC subpopulations of NGC 4365 are 
visible as distinct sequences at g — z ~ 0.9 and ~ 1.3 respec- 
tively. At radii smaller than 2 arcmin and possibly between 
2.5 and 3.5 arcmin there is evidence of a third subpopulation 
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galactocentric radius (arcminutes) 

Figure 15. GC density witii colour and galactocentric radius for the S-Cam data. The colour scale indicates object number per square 
arcmin per colour bin ranging from black at low density to red at high density. The plot is constrained to objects with radii of 0.3 to 
12 arcmin. There is indication of intermediate colour (g' — i' ~ 1.05) GCs at small radii, similar to the ACS data, but it is not clear 
whether this intermediate colour subpopulation is evident at larger radii. 
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bluer colours at larger radii. Forbes et al. (20111 find the 
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Figure 13. Histogram of the azimuthal distribution of GC candi- 
dates. Poisson errors are shown on the histogram and the dotted 
line with best fit ellipticity and position angle is overplotted. See 
text for fitting and parameter details. 



of GCs intermediate in colour. Note the clear overdensity of 
GCs candidates at intermediate to red colours and within 1 
arcmin in both Fig ures|14|and|15| The se data confirm a sim- 
ilar observation by |Larsen et al. ( 2005[ ) where evidence of in- 
termediate colour GCs were found at radii smaller than 1.25 
arcmin. We are able to extend the analysis of the GC colour 
substructure with radius to 5.5 arcmin with HST/ACS imag- 
ing (see Figure 14 1 and to 11.5 arcmin with S-Cam imaging 
(see Figure 15 1. Beyond 3 arcmin in Figure 15 the colour 



distribution seems to be dominated by 2 colour modes but 
there is some indication of GCs at intermediate colours. 

Recently, observers have found gradients in the peak 
colour of the GC subpopulations with galactocentric radius 
(see [Harris] [2009a|b{ [Liu et al.|[20TT] [Forbes et ah] [2011 1. 
Both blue and red subpopulations' peak colour shifts to 



most extreme gradient to date in the galaxy NGC 1407. If 
the red subpopulation of NGC 4365 were similar to that 
of NGC 1407 then the peak colour for the red subpopula- 
tion would change from — i' ~ 1.19 at the galaxy centre 
to — i' ~ 1.06 at 10 arcmin, and the red subpopulation 
would have intermediate colours at large radii. In the bot- 
tom panel of Figure [T4| and in Figure [15| an overdensity of 
GC candidates with intermediate colours is visible at very 
small radii that cannot be accounted for by a gradient in 
the peak colour of the red subpopulation. We note that if 
this intermediate colour overdensity at small galactocentric 
radii is a distinct subpopulation and there is a gradient in 
the red subpopulation of NGC 4365 then statistical tests 
for bimodality or trimodality based on colour distribution 
alone might not be able to distinguish between these two 
subpopulations conclusively. 



6.2 Radial colour gradients 

Several groups have recently found shallow but sigruficant 



gradients in the mean colour of GC subpopulations (see Har 
rls|[2009a|b[ [Liu et al.[[2011[ [Forbes et al.|[2011[ ). The radial 
extent of the Subaru/S-Cam imaging makes this an excellent 
data set to determine whether there are gradients present in 
the GC subpopulations of NGC 4365. The presence of an in- 
termediate colour subpopulation will influence the gradient 
determined and therefore we consider the case of a bimodal 
and a trimodal distribution separately. The heteroscedastic 
KMM algorithm was run on a rolling sample of GC candi- 
dates and the peak value determined for both bimodal and 
trimodal instances for the S-Cam data set. The rolling sam- 
ple consisted of 360 GC candidates input to KMM at each 
step and between each step the 120 candidates with lowest 
galactocentric radius were replaced with the 120 candidates 
closest to, but at least as far as, the furthest candidate in 
the previous sample. A linear relation was fitted to the S- 
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Figure 14. GC density with colour and galactocentric radius 
for the ACS data. The colour scale indicates object number per 
square arcmin per colour bin ranging from black at low density to 
red at high density. The plot is constrained to objects with radii 
of 0.02 to 5.8 arcmin. Top; Measured GC density with colour and 
radius is plotted. Centre; The model of red and blue GC den- 
sity calculated using the Gaussian positions and widths from the 
KMM results for the two subpopulations at large radii (r > 3.5 
arcmin). The proportion of GCs in each Gaussian distribution is 
taken from the KMM results for the whole sample and the Sersic 
radial density profile used is calculated in Section 6.1.1. This was 
scaled to have a similar object density to the data at large radii. 
Bottom; The model in the central panel has been subtracted 
from the density distribution in the top panel to show an ap- 
proximate distribution of GCs that are not well modelled by a 
bimodal colour distribution of GCs. This distribution is not a ro- 
bust measure but it does show indications of significant numbers 
of intermediate colour GCs {g — z ^ 1.1), particularly at small 
radii (< 2 arcmin). 



Cam peak colours, to 14.5 arcmin in the bimodal case and 
7 arcmin in the trimodal case, i.e. 

g — i = a + 6i?(arcmin) (5) 

Results are tabulated in Table |4] The linear blue gradient 
values are consistent within errors but the linear red gradi- 
ent is significantly steeper in the trimodal case. The green 
gradient measured in the trimodal case is consistent with 
zero. 

We also fit a line to the logarithm of the radius nor- 
malised by _Re using = 2.1 arcmin, i.e. 

g-i^c + d\og{R/Re) (6) 

Results are tabulated in Table [5] Here both blue and red 
gradients are consistent within errors for the bimodal and 
trimodal cases. The green subpopulation does not have a 

(I20TT1 
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Figure 16. Rolling peak colour values found by the KMM code 
plotted against galactocentric radius. The bimodal case is plot- 
ted in the top panel and the trimodal case in the bottom panel. 
S-Cam (g' — «') rolling peak values are plotted for blue (dashed), 
green (dash-dotted) and red (dotted) modes in the multimodal 
fits. Galaxy g' — i' colours are plotted as black dots with error- 
bars and the solid black lines show best fit gradients in colour-log 
radius space. 



Subpopulation 


Bimodal Case 


Trimodal Case 




(dex per arcmin) 


(dex per arcmin) 


Blue 


-0.0037 ±0.0002 


-0.004 ± 0.001 


Green 




0.001 ± 0.012 


Red 


-0.0047 ±0.0005 


-0.009 ± 0.002 



Table 4. GC colour gradients for the linear fit to colour with 
radius in arcmin (Equation 5) for S-Cam data. 

tral ACS pointing in the bimodal case, finding a blue slope 
of -0.056 ± 0.024 and a red slope of -0.033 ± 0.021 dex per 
dex. The radial range of their measurements do not over- 
lap with ours and we see that the radial colour gradient 
of the blue subpopulation remains constant when moving to 
larger galactocentric radii while the radial colour gradient of 
the red subpopulation steepens significantly with increasing 
galactocentric radius. We use the |Lee, Park fc Hwang| ( |2010[ ) 
empirical relationship between GC colour and metallicity to 
give a metallicity gradient, i.e. 

A[Fe/H] = 3.48A(g - i) (7) 

The gradient of the mean metallicity with galactocentric 



Subpopulation 


Bimodal Case 


Trimodal Case 




(dex per dex) 


(dex per dex) 


Blue 


-0.055 ±0.004 


-0.04 ± 0.01 


Green 




-0.01 ± 0.10 


Red 


-0.062 ±0.008 


-0.07 ± 0.02 



measurable gradient in this case either. Liu et al 



measure colour gradients for NGC 4365 GCs using the cen- 



Table 5. GC colour gradients for the logarithmic fit to colour 
with normalised radius (Equation 6) for S-Cam data. 
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Subpopulation 


Bimodal 


Case 


Trimodal Case 




(dox per 


dex) 


(dex per dex) 


Blue 


-0.19 ± 


0.01 


-0.13 ±0.03 


Red 


-0.22 ± 


0.03 


-0.26 ±0.06 



Table 6. Metallicity gradients for the logarithmic fit to colour 
with normalised radius (Equations 6 & 7) for S-Cam data. 



radius is shallower in the trimodal case for the blue subpop- 
ulation but for the red subpopulation the gradients in the 
bimodal and trimodal cases are consistent. We note that the 
rolling peak colour found by KMM is only consistent with 



the galaxy colour in the trimodal case (see Figure 161. 

We conclude that the NGC 4365 GC subpopulations 
show a shallow but significant colour gradient with galacto- 
centric radius, regardless of whether the GC system is in- 



terpreted as being bimodal or trimodal. Forbes et al. (2011 1 
tabulate metallicity-radius gradients in units of dex per dex 
for several galaxies on which this measurement has been 
done. The gradients measured for NGC 4365 are steeper 
than the mean, but not the steepest measured. Results vary 
from —0.10 to —0.38 for the blue subpopulation and from 
—0.10 to —0.43 for the red subpopulation in [Forbes et al.| 
( |2011[ ) Table 1 . These gradients explain why statistical anal- 
yses of the colour distribution (collapsed in radius) of NGC 
4365's GC system give inconclusive results. 



6.3 Magnitude colour trends 

In the determination of the magnitude colour trends we only 
investigate the blue subpopulation of the GC system, we 



search for a 'blue tilt' ( jStrader et al.||2006 l. We do this in 
order to decide whether the blue tilt could be the cause of 
the intermediate colour overdensity seen in Figures [4) |14| and 
|15| The bimodal KMM mean colour values at various magni- 
tudes were determined by dividing the GC sample into mag- 
nitude subsets with equal numbers. The mean z magnitude 
of GC candidates in each subset is plotted against the blue 
heteroscedastic Gaussian mean value found by KMM in Fig- 
ure 17 We did this analysis for the central HST/ACS point- 
ing as well as for the full set of eight HST/ACS pointings and 
found that both cases show a significant gradient but that 
the slope is steeper in the central region (—0.0477 ± 0.0018) 
than for all pointings combined (—0.024 ± 0.013). Mieske 



et al. ( 2006 1 use the central ACS pointing and combine GCs 
for all galaxies in the range —21.7 < Mb < —21 to deter- 
mine a slope of — 0.037±0.004 using the same KMM method. 
In Figure[l7]we see that this lies very close to our fitted line. 
They also find a steepening slope when GCs closer the the 
galaxy centre are used. 

In Figure |4j indication of a clump of intermediate 
colour/green objects is seen around g — z ~ 1.1 and 22 > 
2 > 20 and in the bottom panel of Figure [l4] a clear over- 
density is seen inside 2 arcmin and around g — z ~ 1.1. The 
colour magnitude slope for the blue GC candidates is not 
steep enough to cause either of these overdensities. 



6.4 Quantifying the GC subpopulations 

We compare NGC 4365's GC colour distribution with that 
of other Virgo cluster galaxies, using the ACS Virgo Clus- 
ter Survey (VCS) photometry ( [Peng et al.1 '2006'). They fit 
a homoscedastic bimodal distribution using Kaye's Mixture 
Model (KMM) algorithm and obtained mean g — z colour 
values of 0.98 and 1.36 for the blue and red subpopulations 
of NGC 4365 respectively, as well as a width of 0.15 for 
both subpopulations. They also fit linear relations between 
mean g — z colour and absolute galaxy magnitude, Mb, for 
both GC subpopulations, assuming all Virgo Cluster galax- 
ies have the same distance. The relations predict mean g — z 
colour values of 0.98 ± 0.06 and 1.40 ± 0.08 for the blue 
and red subpopulations when the conventional distance for 
NGC 4365 (23.1 Mpc) is used. They concluded that aside 
from having a large number of red GCs, NGC 4365 has an 
unremarkable colour distribution compared to their sample 



( Peng et aL[2006 K 

Fitting a bimodal homoscedastic distribution to NGC 
4365's GCs, with KMM, using all eight available ACS point- 
ings we find similar values to the ones found by ACS VCS 
(i.e. g ~ z — 0.94 and 1.30 for the blue and red peak colours 
and a width of 0.12). However, on visual inspection the fit 
to the GC distribution is skewed to redder colours when 
compared to the actual colour distribution. We also fit both 
bimodal and trimodal heteroscedastic distributions to the 
GC colour distributions, obtaining mean blue and red val- 
ues ol g — z = 0.89 and 1.25 in the bimodal case, and 0.88, 
1.12 and 1.34 in the trimodal case. In neither case does the 
mean blue colour lie close to the g — z vs. Mb relation (val- 
ues are 1.4a and 1.6a away from the relation) but only in 
the trimodal case does the red mean colour lie close to the 
Peng et al. ( 2006 ) g — z v.s. Mb relation, where the value 
is 0.7cr away from the relation instead of 1.8cr away. This 
is a motivation to describe the GC system of NGC 4365 in 
terms of three rather than two subpopulations. 

From this point forward we assume that the GC system 
of NGC 4365 contains three subpopulations, a blue, green 
(intermediate colour) and red subpopulation. There are ra- 
dial gradients in colour and a blue tilt but in neither case 
are they enough to explain the strong overdensity seen at 
intermediate colours in both Figures [l4] and [15] In addition, 
there is better agreement with the colour - galaxy luminos- 



ity relations ofiPeng et al. (20061 to the other Virgo cluster 



galaxies and a near perfect match between galaxy light and 
mean red subpopulation colours if we assume that NGC 4365 
has three subpopulations. 

In the trimodal case, the green GCs would likely 'con- 
taminate' the blue and red subpopulations to all except the 
bluest and reddest colours. A simple split based on colour 
would leave a high percentage of green GCs in both the 
blue and red subpopulations and visa versa. Therefore, we 
formulate a probability that a GC belongs to either a blue, 
green or red subpopulation based on its colour and galacto- 
centric radius. Note that red (metal rich) GCs are generally 
found to be more centrally concentrated than blue (metal 
poor) GCs ( Brodie fc Strader^2006) thus GC galactocentric 
radius is an important factor in the probability calculation. 
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Figure 17. The colour-magnitude relation for blue GC candidates found in the HST/ACS photometry. The left panel shows GC 
candidates found in the central ACS pointing and the right panel shows candidates found in all eight ACS pointings. The KMM mean 
value for the blue GCs is overplotted with blue points. The one sigma range for the linear fit to those values is overplotted with 
a blue shaded region on both panels. On the left panel the gradient the slope found by [Mieske et al.| | |2006[ l for galaxies between 
—21.7 < Mb < —21 is overplotted with a dashed line. 



Table 7. Gaussian values for the blue and red GC distributions as 
determined from an Epanechnikov smoothing kernel of the colour 
distribution. 



green GCs would not be present in significant numbers at 





ACS (g-z) 


S-Cam (g'-i') 




peak width 


peak width 


Blue 


0.89 0.07 


0.80 0.07 


Red 


1.32 0.12 


1.13 0.10 



such blue colours). See Figure 18 for the S-Cam colour distri- 
bution and kernel smoothing. The peak position and width 
of the red distribution were calculated similarly and all val- 
ues are shown in Table [7| As seen in the bottom right panels 
of Figures [19] and [20] where GCs at all radii are included for 
each sample, the Gaussian distributions for blue and red 
subpopulations determined in this way are a reasonable fit 
to the extreme blue and red ends of the GC colour distribu- 
tion. 



We assume that each of the three subpopulations are fairly 
well described by a Gaussian in colour (without skewness 
or kurtosis) and that the peak and width of the Gaussian 
is independent of galactocentric radius (i.e. no radial colour 
gradient) but make no further assumptions about the ra- 
dial distribution of any of the three subpopulations. Given 
evidence of radial colour gradients we note that this assump- 
tion is strictly inaccurate. We make this assumption because 
colour gradients are shallow compared to the width of the 
GC colour distributions and the calculation of gradient val- 
ues contain other significant assumptions about the radial 
distribution of the subpopulations that we wish to avoid. 
The following paragraphs describe how we assign a proba- 
bility that a GC is blue, green or red (henceforth blue-green- 
red probability). 

The peak of the blue Gaussian distribution was de- 
termined from the bluest peak of an Epanechnikov kernel 
smoothing of the colour distribution ( Silverman]] 1986 1 and 
the width of the blue Gaussian distribution was determined 
from the colour at which 68 per cent of the GCs bluewards 
of the peak were included in the distribution (assuming that 



To assess the number of blue, green and red GCs as a 
function of colour and radius, the GCs in each sample (ACS 
or S-Cam) were split into radial bins of roughly equal ob- 
ject number and for each radial bin blue and red Gaussian 
distributions were fitted to the normalised colour histogram. 
Peak and width values were held constant across radial bins, 
as determined from the whole sample (see Table ]7|, and the 
normalisations of the distributions were fitted separately for 
each radial bin. This was a minimization fit. The blue 
and red numbers were calculated by adding the GCs in the 
extreme colour regions (where no significant green 'contam- 
ination' is expected) to the percentage of the total objects 
expected to belong to either the blue or red subpopulations 
(obtained by integrating the blue and red Gaussian distri- 
butions over intermediate colours). The green GC numbers 
were defined by subtracting both the blue and red subpop- 
ulation numbers from the total number in the distribution. 
These values were used to determine a radial surface density 
profile for each subpopulation (discussed further in Section 
6.5.1), which is used to calculate the probability that a GC 



© 0000 RAS, MNRAS 000, 000-000 



Wide-field imaging of NGC 4365's globular cluster system: The third subpopulation revisited 15 




0.92'<r<1.5' 

H 1 1 h 




H h 




4.2'<r<5.8' 



all radii 



0.8 0.9 1.0 1.1 1.2 1.3 1.4 1.5 

g-z 



0.8 0.9 1.0 1.1 1.2 1.3 1.4 1.5 

g-z 



Figure 19. The normalised distribution of ACS GCs with colour in the seven radial bins (shown in arcminutes on each plot) and at 
all radii (lower right panel). Histogram heights are shown as black circles with errorbars representing Poissonian errors. Fitted Gaussian 
distributions are plotted in blue and red for the respective subpopulations and the excess GCs (calculated only in the region between 
the blue and red peaks) are plotted as green squares. The peak and width parameters for the Gaussian distributions fitted in each case 
are constant: fii, = 0.89, ui, = 0.07, fir = 1.32 and ar = 0.12, with the normalisation allowed to vary. 
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Figure 20. The normalised distribution of S-Cam GCs with colour in seven radial bins (shown in arcminutes on each subplot) and at 
all radii (lower right panel). Histogram heights are shown as black circles with errorbars representing Poissonian errors. Fitted Gaussian 
distributions are plotted in blue and red for the respective subpopulations and the excess GCs (calculated only in the region between 
the blue and red peaks) are plotted as green squares. The peak and width parameters for the Gaussian distributions fitted in each case 
are constant: fii, = 0.80, ui, = 0.07, fir = 1.13 and ar = 0.10, with the normalisation allowed to vary. 
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Figure 21. GC radial surface density for blue (squares), green 
(triangles) and red (circles) subpopulations incorporating ACS 
(unfilled points) and S-Cam (filled points) GCs. The surface den- 
sity for both the blue and red GC subpopulations are well fit by 
a Sersic profile added to a background term (shown as blue and 
red solid lines). The surface density for the green subpopulation 
is best fit by a power law plus a background value (shown as a 
green solid line). The red dashed line is a scaled version of the 
Sersic profile fit to the red subpopulation that is a reasonable fit 
to all except the inner two green surface density points. See text 
for further details. 



Figure 18. Colour distribution of GC candidates at all radii 
for S-Cam (top) and ACS (bottom) brighter than the turnover 
magnitude {z = 23 A and i' = 23.6). The Epanechnikov kernel 
smoothing of the colour distribution is plotted with a dotted line, 
the smoothing width used is the average colour error (0.04 for 
both ACS and S-Cam). At least three peaks can be seen. 







n 




bg 




(arcmin~^) 




(arcmin) 


(arcmin"^) 


Blue 


1.67 ± 0.29 


1.36 ±0.19 


7.30 ± 0.68 


0.25 


Green 


1.71 ± 0.17 


2.02 


3.17 


0.21 ± 0.03 


Red 


4.91 ± 0.59 


2.02 ±0.25 


3.17 ±0.19 


0.50 



at any radius or colour belongs to the blue, green or red 
subpopulations. 

The process used to determine the blue-green-red prob- 
ability of a GC starts by determining the relative number of 
GCs in each subpopulation at the galactocentric radius of 
the GC. This is done by comparing the values of the radial 
surface density profile for each subpopulation. The relative 
numbers of GCs in each population are used to scale the 
normalisations of the three Gaussian distributions and com- 
paring the relative values of all three Gaussian distributions 
at the colour of the GC we calculate the probability of an 
object being blue, green, or red. 

6.5 Characterising the GC system subpopulations 

6.5.1 Surface density 

As described in Section 6.4 we counted the blue, green and 
red GCs in each radial bin assuming that the red and blue 
GC colour distributions are well fit by Gaussians. In Figure 
|21| the radial surface density data for all three subpopu- 
lations are shown with the best fit surface density model 
in each case. The red and blue subpopulation distributions 
are both well fit by a Sersic profile added to a background 



Table 8. Sersic profiles fits to the radial surface density of the 
blue, green and red GC subpopulations. The green subpopulation 
n and Re values have been fixed to the red subpopulation values. 

term. The values we found are recorded in Table d The 
background terms were estimated from the outermost radial 
surface density value. The effective radius, or half number 
radius, of the blue subpopulation (7.30 arcmin) is more than 
twice as large as the effective radius of the red subpopulation 
(3.17 arcmin) and the Sersic n values are consistent within 
combined errors. The green subpopulation can be fit by the 
red subpopulation Sersic profile if the two innermost density 
points are omitted (points interior to 1.5 arcmin) and we fit 
for Pe and the bg (see Table [s] for the fitted values). 

We also fit a two parameter power law plus a back- 
ground term 

p{R) = poR" + bg (8) 

to the entire radial range of the green subpopulation, find- 
ing a = —1.71 ± 0.13, po = 12.3 ± 1.3 arcmin"^ and 
bg = 0.10 ± 0.05 arcmin"'^. Comparing the slope of a power 
law fit to the blue (-1.13 ± 0.06) and red (-1.60 ± 0.04) 
distributions over a radial range from 0.5 to 11 arcmin (4 to 
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Figure 22. The percentage contribution of eacii subpopulation 
to the total number of GC candidates plotted against galacto- 
centric radius. Curves were calculated using the fitted surface 
density profiles in Figure |21| Both the data points and plotted 
profiles have had the respective background values subtracted. 
Beyond 12 arcmin the total population numbers are low and the 
measurement of the percentage contribution of any one subpop- 
ulation to the total number becomes unstable. We see that the 
red subpopulation dominates at small galactocentric radii and the 
blue dominates at large galactocentric radii while the green sub- 
population is only significant at very small galactocentric radii. 



70 kpc) we can quantitatively verify that the green and red 
GC subpopulations are more centrally concentrated than the 
blue GC subpopulation and that the slope of the green and 
red subpopulations are consistent with each other. 

In Figure [22] the percentage of the total GC candidate 
number that is attributed to each subpopulation is plotted 
as a function of galactocentric radius. The data points have 
been background subtracted and the curves derived from 
the profile fits in Figure [21] are plotted without the addition 
of background terms. In the case of the green subpopulation 
the curve shown is a combination of the power law fit in 
the inner parts {R < 2.25 arcmin) and the Sersic profile be- 
yond that. We do this to adequately display the behaviour 
of the intermediate colour surface density both in the inner 
(power law behaviour) and outer (Sersic behaviour) parts, 
noting that we plot lines to guide the eye and not for analy- 



sis. In agreement with previous work (see review by Brodie & 
[Strader^2006J the red subpopulation dominates the GC sys- 
tem in the inner parts and the blue subpopulation dominates 
in the outer parts. We also see that the green subpopulation 
only contributes significantly to the GC system at very small 
galactocentric radii {r < 2 arcmin). This is in agreement 
with the qualitative analysis of the colour-galactocentric ra- 
dius properties of the GC system earlier in this work (see 



Section 6.1.2) and with Larsen et al. (20051 



6.5.2 GC half light sizes 

Here we measure subpopulation sizes (median half light 
radii) as well as trends in candidate size with galactocentric 
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Figure 23. The half light radius (r^) of HST/ACS GC can- 
didates against galactocentric radius for objects brighter than 
z = 23.4. The high probability blue, green and red GC sub- 
population candidates are plotted with blue squares, green tri- 
angles and red circles respectively. There is a trend of increasing 
GC candidate size with increasing galactocentric radius for all 
GCs regardless of subpopulation colour. A power law fitted to 
the whole sample of GC candidates brighter than z = 23.4 has a 
slope of 0.49 it 0.04 dex per dex. Overplotted with blue dashed, 
green dash-dotted and red dotted lines are the power law fits for 
size normalisation (with slope set by whole GC sample) of each 
high probability subpopulation sample. 



radius from the ACS photometry. As described in Section 6.4 
each GC candidate is assigned a probability of belonging to 
each subpopulation. For example a candidate 1.28 arcmin 
from the galaxy centre with g — i = 0.98 has probabilities: 
Pbiuc = 18 per cent, Pgrccn = 79 per cent and Prod = 3 per 
cent. Defining the blue and red GCs with a colour probabil- 
ity larger than 95 per cent and green GCs with a probability 
larger than 80 per cent there are 315, 426 and 491 GC can- 
didates in each of the blue, green and red high probability 
subsamples. There are no green GC candidates with a prob- 
ability larger than 95 per cent so we compare medium prob- 
ability green candidates (80-95 per cent) with high proba- 
bility blue and red candidates. The aforementioned example 
is not defined as part of either blue, green or red high prob- 
ability subsamples. 

We determine the median half light radius of the GCs; 

?-h,biuc = 0.037 ±o:o!]i arcsec 

n OOT I 0.001 

ni.grocn = 0.027 ±0.001 arcsec 

n nnr: i 0.001 

»"h,rcd = 0.025 ±0.001 arcsec 

These sizes were measured on GCs brighter than z = 23.4 
(109, 223 and 269 blue, green and red high probability 



subsample GCs) using ISHAPE ( Larsen 1999 1 where the 



method is described in Strader et al. (20061. The median 



sizes correspond to 4.1±o:2, 3.0±o:i and 2.8±o:i parsec re- 
spectively. Blue GCs have a significantly larger median size 
than either green or red GCs. However, the sizes of the green 
GCs are only slightly larger than the red GCs. This is pos- 
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Figure 24. The normalised mass function for GC candidates 
brighter than z = 23.4 from the ACS catalogue. All GCs (black) 
as well as the high probability subpopulation candidates (plot- 
ted in blue dashed, green dash-dotted and red dotted lines re- 
spectively). Each subpopulation has a significantly different mass 
distribution to the others. 



sibly due to distinct characteristic sizes for each subpopu- 
lation but could also be explained by a continuous trend in 
GC size with GC colour (e.g. |Jordan|2004 |. 

As well as a possible trend of GC size with colour we 
also see a trend of GC size with galactocentric radius. In Fig- 
ure [23] a trend of increasing GC candidate size with increas- 
ing galactocentric radius is visible. Fitting a two parameter 
power law to the GC candidates brighter than z = 23.4 we 
find 



r-h(pc) = [1.00±0.04]i?(kpc 



[0.49±0.04) 



All three subpopulations show a clear size increase with 
distance from the galaxy centre with a power law slope 
of 0.49 ± 0.04 dex per dex. This measurement extends to 
4.5i?e, one of only four similar measurements extending be- 
yond ~ 2i?e. |Harris| ( [2009a) measured a slope of 0.11 on a 
composite sample of six massive galaxies, [Gomez fc Wood-] 
ley (20071 measured a slope of 0.05 ± 0.05 for metal poor 



GCs and 0.26 ± 0.06 for metal rich GCs in NGC 5128 and 



Spitler et al. (20061 measured a slope of 0.19 ±0.03 in NGC 



4594. The slope we measure is steeper than any of the pre- 
vious values and much closer to the value measured for our 



own Galaxy (0.36 ± 0.07 for metal rich GCs, see Gomez & 
Woodley]]2007 1 . We conclude that there is more variability. 



between galaxies, in the relationship of GC size with galac- 
tocentric radius than previously found. Shown in Figure [23] 
are the GCs in all three high probability subpopulation sam- 
ples and power law fits. We set the slope of the fits to 0.49 
and found normalisations of 1.096±0.042, 0.994±0.059 and 
0.944±0.044 for the blue, green and red samples respectively. 
This results agrees with analysis of GC subpopulation me- 
dian sizes. 



6.5.3 Mass function 

We compare the mass function of the high probability sub- 
population GC candidates (as defined in Section 6.5.2). The 
mass for each GC candidate was derived using the colour 
dependent values for the z filter mass-to-light (M/L^) ratio 



in Table 5 of Jordan et al. (20071. They used PEGASE 2.0 



with a modified Salpeter initial mass function from |Ken-| 
(19831 and an age of 13 Gyr for all GCs. We used 



nicutt 



magnitudes brighter than z — 23.4 to calculate the mass 
functions (plotted in Figure [24|. We ran the Kolmogorov- 
Smirnov (KS) test on each pair with a common mass range 
of 2 X lO^M© < M < IO'^Mq. The KS test gives a probabil- 
ity of 98.3 per cent that the blue and red subpopulations are 
drawn from different distributions and gives a probability of 
more than 99.9 (99.4) per cent that the green subpopulation 
is drawn from a different distribution than the blue (red) 
subpopulation. The discrepancy between the three subpop- 
ulations is greatest around 1 x IO^Mq in the cumulative 
mass distribution the KS test uses for its analysis. 

Analysis of the mass function of GC subpopulations 
indicates that the green subpopulation is distinct from 
both blue and red subpopulations, containing a significantly 
larger percentage of objects between IO^Mq and lO^M© 
than either blue or red subpopulations, as can be clearly 
seen in Figure |24| This overdensity of green subpopulation 
objects can also be seen in Figure [4] around g — z ~ 1.1 and 
22 > z > 19. It is possible that UCDs or dE nuclei contam- 
inate the bright, mid-colour range ( [Strader et al. 20061 to 
produce this feature. 



6.5.4 Azimuthal properties 

In Figure |25] we show the azimuthal distribution of high 
probability colour subpopulation GC candidates (as defined 
in Section 6.5.2) for S-Cam and ACS. We restrict the sam- 
ples to radii where we have complete azimuthal coverage so 
as not to bias the samples in any direction. The position an- 
gle data is folded at 180° to improve number statistics. The 
bottom panel of Figure[25]shows that there is no position an- 
gle signal within errors for any of the three subpopulations in 
the ACS sample. The S-Cam azimuthal distribution in the 
top panel does show some structure; the high probability 
blue GC candidates are more likely to have position angles 
between 0° and 90°, which is consistent with the blue GC 
subpopulation being elongated along the major axis of NGC 
4365 (the galaxy has a position angle of ~ 42°, see Section 
5). The red distribution shows an almost fiat azimuthal dis- 
tribution (indicating an almost circular GC system) and the 
green distribution shows a very shallow sinusoidal distribu- 
tion with a peak between 0° and 90°. It is clear that the blue 
GC system is very elongated along the position angle of the 
galaxy and likely that both green and red subpopulations 
also have a similar position angle to the galaxy starlight. 
We obtain estimates of the ellipticity of the blue, green and 
red subpopulations by fixing the position angle (to 42°) as 
well as the power law exponent and background value as de- 
scribed in Section 6.1.1, Equation 4 and following the same 
fitting procedure. We find that the estimated ellipticity for 
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Figure 25. Histograms of the azimuthal distributions of liigh 
probability blue (dashed), green (dash-dotted) and red (dotted) 
GC subpopulations. The top panel shows the position angles for 
the S-Cam GC candidates out to 11 arcmin from the galaxy cen- 
tre and the bottom panel shows the ACS GC candidates out to 
3.38 arcmin from the galaxy centre. The position angle of the 
galaxy is ~ 42° and both blue and green subpopulations show 
GC overdensity between 0° and 90° in the S-Cam catalogue. 
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Figure 26. The radial surface brightness profile of the galaxy 
light from the S-Cam i' filter compared with the surface density 
profile of the GC subpopulations. The GC subpopulations are 
plotted as blue squares, green triangles and red circles. The grey 
points show the surface brightness measurements and the black 
line shows the Sersic profile fit to the galaxy light (see Section 
5.1). All parameters are plotted against equivalent galactocentric 
radius, calculated as the geometric mean of the semi-major and 
semi-minor radii. The green and red GC surface density profiles 
are similar to the galaxy light surface brightness profile. 



The agreement in ellipticity between the galaxy light 
and the red GC subpopulation is consistent with literature 
findings that red GC subpopulations generally follow the 
properties of galaxy field stars more closely than blue GC 
subpopulations. The blue GC subpopulation dominates in 
the outer regions of the galaxy and also dominates the mea- 
surement of the ellipticity of the GC system. It can be seen 
qualitatively in Figure [12] and quantitatively from the re- 
sults in Table [9] that the blue subpopulation is significantly 
more elliptical than the galaxy field stars. 



Galaxy Light Blue Green Red 

0.25 ±0.03 0.66 ±0.06 0.55 ± 0.07 0.16 ±0.25 



Table 9. Estimated ellipticity for blue, green and red GC sub- 
populations. The position angle of each GC subpopulation is fixed 
to be the same as that of the galaxy light (42°). 

the blue and green subpopulations are the same within er- 
rors and also consistent with the ellipticity measured for the 
whole GC sample in Section 6.1.1. The estimated ellipticity 
of the red subpopulation is much smaller than the blue and 
green values. It is consistent with zero and also the ellipticity 
measured for the galaxy light (see Table [9]|. 

We also split each subpopulation into an inner and outer 
radial bin with equal numbers in each and refit for ellipticity 
but did not detect radial variation in ellipticity for any of 
the three subpopulations. 



6.6 Comparison with galaxy surface brightness 

The galaxy surface brightness profile from the i' filter pho- 
tometry is compared with the GC surface density profiles in 
Figure |26[ We compare the galaxy surface brightness pro- 
file with 2.5 times the logarithm of the GC surface density 
profile and arbitrarily scaled for ease of viewing, i.e. 

2.51og[P(i?) - bg] (9) 

where P{R) is the GC radial surface density profile and bg 
is the determined background value. Visually, the shape of 
the blue GC surface density profile is very different to that 
of the galaxy surface brightness whereas both red and green 
GC surface density profiles have shapes similar to the galaxy 
light. Both red and green density profiles are slightly steeper 
than the galaxy light profile beyond ~ 5 arcmin. 
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7 DISCUSSION 

In previous work on NGC 4365, evidence of an additional 
subpopulation of GCs at intermediate (green) colours was 
based on small numbers and restricted to the central regions 



(|Puzia et al.| 2002 


Larsen et al. 


2003 


Hempel & Kissler- 


Patig||2004| IBrodie et al.||2005| |Kundu et al.||2005| |Larsen 


et al.|2005l. The lar 


ge sample of GCs in this work, covering 



an extended radial range, allows us to revisit these claims 
in the optical wavelength range. 

We applied several independent statistical tests to the 
colour distribution, i.e. Chi-Squared minimization (x^), the 
Kolmogorov-Smirnov (KS) test and Kaye's Mixture Model 
(KMM) algorithm. When the entire radial range of our imag- 
ing is used our statistical tests cannot conclusively rule out a 
bimodal interpretation, partly because GC subpopulations 
do not have strictly Gaussian colour distributions. 

When examining GC subpopulation colours with ra- 
dius we find significant negative colour/metallicity gradi- 
ents in both the bimodal and trimodal cases. Due to its 
radial colour gradient, the red subpopulation has somewhat 
greener colours at large radii. This explains why statistical 
tests showed inconclusive results when the full radial ex- 
tent of the imaging was used. However, this cannot explain 
the central green overdensity that motivates splitting NGC 
4365's GCs into three subpopulations. 

We compared the photometrically observable proper- 
ties of blue, green and red GCs (measured here to contain 
43, 17 and 40 per cent of total GCs respectively). The high 
probability blue, green and red samples have colour ranges 
oi g - i < 0.82, 0.89 < g - i < 1.0 and g - i > 1.1 where 
definitions are based on kernel smoothing of the GC colour 
distribution. We find evidence for different properties be- 
tween the subpopulations. The green GCs are intermediate 
in size to the blue and red subpopulations (clearly smaller 
than the blue GCs and marginally larger than the red ones) 
and have a significantly difi'erent mass function to both blue 
and red subpopulations. The size differences could be due 
to a continuous trend with colour rather than specific sub- 
population properties (see |Jordan||2004| ) but the differences 
in mass functions suggest that the three subpopulations are 
indeed distinct. Green GCs also have a significantly steeper 
surface density profile than red GCs within 2 arcmin (13.4 
kpc) from the galaxy centre. However, the green and red sub- 
populations can be described by the same radial surface den- 
sity profile, with different normalisations, beyond 2 arcmin 
(13.4 kpc). The ellipticity of the green GC subpopulation is 
similar to the blue GC subpopulation and significantly more 
elongated than that of the red. We do not detect any varia- 
tion in ellipticity of the subpopulations with galactocentric 
radius to 11 arcmin (74 kpc). These results lend support to 
the idea that blue, green and red GCs form three distinct 
subpopulations. 

With the same three subpopulation division we find 
a number of similarities between the properties of the red 
subpopulation and the galaxy light of NGC 4365. The g — i 
radial colour gradient of the galaxy light and mean colour 
of the red GC subpopulation are virtually identical (see 
also ISpitlerl |2010| iForbes & Fortel 120011). The eUipticity 



of the red subpopulation is consistent with that of the 
galaxy starlight and the surface density profile of the red 
subpopulation is closer to the starlight surface brightness 
profile than either blue or green subpopulations. These 



results confirm similar findings in the literature (Brodie & 
Strader|2006 1 and support GC formation scenarios in which 



the red subpopulation is closely linked to the formation of 
the bulk of the galaxy field stars. 

To understand the physical implications of optical 
colour trimodality we consider the relationship between op- 
tical colours, metallicities and ages. There is considerable 
debate about whether the colour distribution of GC systems 
is a good proxy for the metallicity distribution, even if uni- 
form old ages for GCs are assumed (see discussions on the 
non-linearity of colour-metallicity tran sformations by | Yoon| 



|et al.|2006[|Cantiello fc Blakeslee|2007[ [Hempel et al.|2007 l 
It is also possible that one or more of the GC subpopulations 
have significantly younger ages caused by a merger (see 



Hi- 



bbard fc Mihos|[T995l [Schweizer fc Seitzer|[T998l [Montuori 
et al.pOlO l) though there is little spectroscopic evidence of 



young ages for NGC 4365 GCs | |Brodie et alT][2005 K 

The green subpopulation might be a consequence of a 
unique evolutionary history of NGC 4365. The kinematically 
distinct core (KDC) in the stellar light might be a signal of 
this. This property is only found in galaxies the SAURON 
team ( de Zeeuw et al.||2002 | classify as slow rotators in the 
inner regions ( [Emsellem et al.|[2007l ). The KDC is confined 
within 5 arcsec (0.56 kpc) of the galaxy centre and the stellar 
population inside and outside the KDC are indistinguishable 
in age and metallicity. Inside 5 arcsec the stars rotate around 
the minor axis and outside 5 arcsec they rotate around the 
major axis ( Davies et al.|2001 l. Davies et al. (20011 suggest 
that the KDC might have been formed in the merger of "gas- 
rich fragments at high redshift" and state that they find no 
evidence of recent star formation in NGC 4365. The effective 
radius of the galaxy (~ 2 arcmin or ~ 13.4 kpc) and the 
radial extent of the green subpopulation (~ 2 arcmin) are 
much larger than the KDC. It is possible that the formation 
of the green GC subpopulation is linked to the formation 
of the KDC but since they are not spatially correlated, and 
other elliptical galaxies with KDCs show no clear indication 
of a green subpopulation, the connection is not clear. 

The core must have been formed very early on because 



Davies et al. (20011 find the same > 12 Gyr luminosity- 
weighted age for the KDC and the rest of NGC 4365, with 
no sharp changes in metallicity across the boundary of the 
KDC. If the formation of the green subpopulation is linked 
then we would expect the green subpopulation to also be 
very old. [Brodie et air| ( |2005[ ) found evidence that the green 
subpopulation was indeed old but this analysis was based 
on a small number of GCs. We might also expect to find 
rotation for the green subpopulation to be offset by ~ 90° 
compared to the other GC subpopulations and the bulk of 
the starlight, if it is associated with the formation of the 
KDC. 

Another possible signature of the unique evolutionary 
history could be the significant misalignment between the 
galaxy's kinematic and photometric major axes. While 



© 0000 RAS, MNRAS 000, 000-000 



22 Christina Blom, Lee R. Spitler and Duncan A. Forbes 



KDCs are relatively common in elliptical galaxies, the 
minor axis rotation of the bulk of NGC 4365's starlight 
is relatively uncommon. Van den Bosch et al. (2008) use 
triaxial orbit based models to explain NGC 4365's apparent 
minor axis rotation and KDC. The combination of axisym- 
metry in the inner parts with triaxiality in the outer parts 



that van den Bosch et al. ( 2008 1 describe is consistent with 



the age and metallicity of the KDC but bears no obvious 
relation to the formation of a green GC subpopulation 
in NGC 4365. Alternatively, [Hoffman et al.| ( |2010[ ) show 
that a ~ 2.5 Gyr old major merger remnant with 15 — 20 
per cent progenitor gas fraction has a kinematic signature 
remarkably similar to NGC 4365. The gas rich major 



merger that Hoffman et al. (20101 describe is a tempting 



explanation for the presence of green GCs because all NGC 
4365's anomalies would have one explanation. Given that 
the green GC subpopulation is likely old ( Bro die et al.|2005 1 
and there is no evidence for young stars in the KDC ( Davies 
et al.|200T I this would have to be a very early major merger. 



The current galaxy and GC system formation scenar- 
ios that describe mechanisms for producing bimodality in 
GC systems could also explain GC system trimodality for a 
subset of galaxies under certain conditions. 

In the multiphase collapse scenario GC systems are 
formed during two phases, the first phase produces metal 
poor GCs, is truncated by one of a variety of possible mech- 
anisms and then later the second phase produces the metal 
rich GCs (Forbes et al. 19971. In this context a GC sys- 



tem that is trimodal in metallicity could simply be the re- 
sult of three formation phases truncated twice (by the same 
or two different mechanisms). The truncation of the metal 
poor formation phase could be caused by a universal epoch 
of reionization but, for truncation to happen twice, at least 
one truncation mechanism has to be related to nearby galax- 
ies or a process internal to the galaxy, like an active galactic 
nucleus. 

The major merger formation scenario ( |Zepf fc Ashman] 



1993), where metal poor GCs are present in galaxies be- 



fore mergers and metal rich GCs are formed during major 
mergers, could be consistent with some galaxies containing 
three GC subpopulations. If a galaxy underwent a very early 
gaseous major merger and later another major merger it may 
have two subpopulations more metal rich than the metal 
poor GC subpopulation. Depending on the age and metal- 
licity difference it could appear trimodal in colour. 



The accretion model for GC system formation (Cote 



et al. 1998 \ found more than two modes in some of their 



simulated GC systems before a third GC subpopulation was 
considered for NGC 4365 from observations. ICote et al.l 



( 1998 1 mention that more than two peaks are present in 



some of their simulations when a very steep luminosity func- 
tion (Schechter function a = —1.8) is used. They find that 
the peak of the metal poor population correlates with the 
slope of the luminosity function and that for very steep 
slopes the metal poor peak is more metal poor. It is con- 
ceivable that the merger histories of some giant elliptical 
galaxies will therefore show the presence of a third, old, in- 
termediate metallicity subpopulation. 



The questions that remain include distinguishing which 
of these formation scenarios best explain the GC system 
of NGC 4365 and determining whether the trimodality in 
the colour distribution of the GC system is reflected in the 



metallicity distribution. The method in Foster et al. (20101 



for determining metallicity from the Calcium triplet indices 
could provide a large sample of metallicity measurements 
over a wide field of view to assess whether the metallicities 
of the three colour subpopulations we define here are dis- 
tinct. The spectroscopic analysis of a large number of GCs 
could also be used to build a picture of the kinematics of 
the GC system of NGC 4365 and its subpopulations. An 
understanding of NGC 4365's GC system kinematics could 
distinguish whether two or three subpopulations were all 
formed in situ, whether the more metal rich populations 
were accreted or whether the subpopulations are the result 
of several major mergers. 



8 SUMMARY AND CONCLUSIONS 

Combining the photometric depth, size information and res- 
olution of HST/ACS data with the spatial extent and three 
filter imaging of the Subaru/S-Cam data we can achieve a 
uniquely detailed, and unmatched spatially extensive, anal- 
ysis of NGC 4365's GC system. The GC system extends 
beyond 134 kpc from the galaxy centre to 9.5-Re. The blue 
GC subpopulation has not yet reached the background level 
at the very edges of our Subaru/S-Cam imaging. We place 
a lower limit on the total number of GCs to be 6450 ± 110. 

We find further evidence to support the existence of a 
distinct third subpopulation at intermediate colours in the 
GC system of NGC 4365. We also find a trend of increasing 
GC size with galactocentric radius and negative gradients in 
the colour/metallicity of both blue and red subpopulations 
with galactocentric radius. The blue subpopulation shows 
evidence for a blue tilt. Comparing the GC system with the 
galaxy light we find that the red subpopulation has a similar 
colour, radial colour gradient, ellipticity and radial surface 
density slope to the galaxy light. 

Based on various measured properties we find it most 
likely that the green subpopulation is distinct from the blue 
and red subpopulations. Consequently, it is very likely that 
NGC 4365 has had a unique evolutionary history causing the 
existence of a third GC subpopulation in this giant elliptical, 
where most similar galaxies only have two. The photomet- 
ric properties of the blue and red GC subpopulations are 
consistent with the properties found for those of other giant 
elliptical GC systems, therefore any formation scenario for 
the green subpopulation must leave the predicted properties 
of the other two subpopulations relatively unchanged. 

Future analysis of spectroscopy already obtained for 
this system will help determine the origin of the green sub- 
population of NGC 4365's GC system. 



9 ACKNOWLEDGEMENTS 

We thank the reviewer, A. Graham, C. Foster, M. Owers, A. 
Romanowsky and A. Chies-Santos for insightful discussions. 



© 0000 RAS, MNRAS 000, 000-000 



Wide-field imaging of NGC 4365's globular cluster system: The third subpopulation revisited 23 



We also thank M. Smith, V. Pota, C. Usher and S. Kartha 
for support during the preparation of this manuscript. Based 
in part on data collected at Subaru Telescope, which is op- 
erated by the National Astronomical Observatory of Japan. 
Based in part on observations made with the NASA/ESA 
Hubble Space Telescope, and obtained from the Hubble 
Legacy Archive, which is a collaboration between the Space 
Telescope Science Institute (STScI/NASA), the Space Tele- 
scope European Coordinating Facility (ST-ECF/ESA) and 
the Canadian Astronomy Data Centre (CADC/NRC/CSA). 



REFERENCES 

Alves-Brito A., Hau G. K. T., Forbes D. A., Spitler L. R., 

Strader J., Brodie ,J. P., Rhode K. L., 2011, ArXiv c-prints 
Arnold J. A., Rornanowsky A. J., Brodie J. P., Chomiuk 

L., Spitler L. R., Strader J., Benson A. J., Forbes D. A., 

2011, ApJ, 736, L26+ 
Ashman K. M., Bird C. M., Zepf S. E., 1994, AJ, 108, 2348 
Ashman K. M., Zepf S. E., 1992, ApJ, 384, 50 
Beasley M. A., Baugh C. M., Forbes D. A., Sharpies R. M., 

Frenk C. S., 2002, MNRAS, 333, 383 
Bekki K., Yahagi H., Nagashima M., Forbes D. A., 2008, 

MNRAS, 387, 1131 
Bender R., Burstein D., Faber S. M., 1992, ApJ, 399, 462 
Blakeslee J. P., Cantiello M., Peng E. W., 2010, ApJ, 710, 

51 

Blakeslee J. P., Jordan A., Mei S., Cote P., Ferrarese L., 

Infante L., Peng E. W., Tonry J. L., West M. J., 2009, 

ApJ, 694, 556 
Brodie J. P., Strader J., 2006, ARA&A, 44, 193 
Brodie J. P., Strader J., Denicolo G., Beasley M. A., Ce- 

narro A. J., Larsen S. S., Kuntschner H., Forbes D. A., 

2005, AJ, 129, 2643 
Burstein D., Davies R. L., Dressier A., Faber S. M., Stone 

R. P. S., Lynden-Bell D., Terlevich R. J., Wegner G., 1987, 

ApJS, 64, 601 
Cantiello M., Blakeslee J. P., 2007, ApJ, 669, 982 
Chen C, Cote P., West A. A., Peng E. W., Ferrarese L., 

2010, ApJS, 191, 1 

Chies-Santos A. L., Larsen S. S., Kuntschner H., Anders P., 
Wehner E. M., Strader J., Brodie J. P., Santos J. F. C, 

2011, A&A, 525, A20+ 

Cote P., Marzke R. O., West M. J., 1998, ApJ, 501, 554 
Davies R. L., Kuntschner H., Emsellem E., Bacon R., Bu- 
reau M., CaroUo C. M., Copin Y., Miller B. W., Monnet 
C, Peletier R. F., Verolme E. K., de Zeeuw P. T., 2001, 
Ap.I, 548, L.33 

de Zeeuw P. T., Bureau M., Emsellem E., Bacon R., Car- 
olio C. M., Copin Y., Davies R. L., Kuntschner H., Miller 
B. W., Monnet C, Peletier R. F., Verolme E. K., 2002, 
MNRAS, 329, 513 

Emsellem E., Cappellari M., Krajnovic D., van de Ven 
G., Bacon R., Bureau M., Davies R. L., de Zeeuw P. T., 
Falcon-Barroso J., Kuntschner H., McDermid R., Peletier 
R. F., Sarzi M., 2007, MNRAS, 379, 401 



Ferrarese L., Cote P., Jordan A., Peng E. W., Blakeslee 
J. P., Piatek S., Mei S., Merritt D., Milosavljevic M., 
Tonry J. L., West M. J., 2006, ApJS, 164, 334 

Forbes D., Spitler L., Strader J., Romanowsky A., Brodie 
J., Foster C, 2011, ArXiv e-prints 

Forbes D. A., 1996, AJ, 112, 954 

Forbes D. A., Brodie J. P., Grillmair C. J., 1997, AJ, 113, 
1652 

Forbes D. A., Forte J. C, 2001, MNRAS, 322, 257 
Forbes D. A., Franx M., lUingworth G. D., Carollo C. M., 

1996, ApJ, 467, 126 
Foster C, Forbes D. A., Proctor R. N., Strader J., Brodie 

J. P., Spitler L. R., 2010, AJ, 139, 1566 
Foster C, Spitler L. R., Romanowsky A. J., Forbes D. A., 

Pota v., Bekki K., Strader J., Proctor R. N., Arnold J. A., 

Brodie J. P., 2011, MNRAS, pp 906—1- 
Gomez M., Woodley K. A., 2007, ApJ, 670, L105 
Goudfrooij P., Hansen L., Jorgensen H. E., Norgaard- 

Nielsen H. U., de Jong T., van den Hoek L. B., 1994, 

A&AS, 104, 179 
Graham A. W., Driver S. P., 2005, PASA, 22, 118 
Harris W. E., 2009a, ApJ, 699, 254 
Harris W. E., 2009b, ApJ, 703, 939 
Harris W. E., van den Bergh S., 1981, AJ, 86, 1627 
Hempel M., Kissler-Patig M., 2004, A&A, 428, 459 
Hempel M., Kissler-Patig M., Puzia T. H., Hilker M., 2007, 

A&A, 463, 493 
Hibbard J. E., Milios J. C, 1995, AJ, 110, 140 
Hoffman L., Cox T. J., Dutta S., Hernquist L., 2010, ApJ, 

723, 818 

Jordan A., 2004, ApJ, 613, L117 

Jordan A., McLaughlin D. E., Cote P., Ferrarese L., Peng 

E. W., Mei S., Villegas D., Merritt D., Tonry J. L., West 

M. J., 2007, ApJS, 171, 101 
Jordan A., Peng E. W., Blakeslee J. P., Cote P., Eyhera- 

mendy S., Ferrarese L., Mei S., Tonry J. L., West M. J., 

2009, ApJS, 180, 54 
Kennicutt Jr. R. C, 1983, ApJ, 272, 54 
Kormendy J., Fisher D. B., Cornell M. E., Bender R., 2009, 

ApJS, 182, 216 
Kundu A., Whitmore B. C, 2001, AJ, 121, 2950 
Kundu A., Zepf S. E., Hempel M., Morton D., Ashman 

K. M., Maccarone T. J., Kissler-Patig M., Puzia T. H., 

Vesperini E., 2005, ApJ, 634, L41 
Larsen S. S., 1999, A&AS, 139, 393 

Larsen S. S., Brodie J. P., Beasley M. A., Forbes D. A., 
Kissler-Patig M., Kuntschner H., Puzia T. H., 2003, ApJ, 
585, 767 

Larsen S. S., Brodie J. P., Huchra J. P., Forbes D. A., 

Grillmair C. J., 2001, AJ, 121, 2974 
Larsen S. S., Brodie J. P., Strader J., 2005, A&A, 443, 413 
Lee M. G., Park H. S., Hwang H. S., 2010, Science, 328, 

334 

Lee M. G., Park H. S., Hwang H. S., Arimoto N., Tamura 

N., Onodera M., 2010, ApJ, 709, 1083 
Liu C, Peng E. W., Jordan A., Ferrarese L., Blakeslee J. P., 

Cote P., Mei S., 2011, ApJ, 728, 116 
McLaughlin D. E., Harris W. E., Hanes D. A., 1994, ApJ, 

422, 486 



© 0000 RAS, MNRAS 000, 000-000 



24 Christina Blom, Lee R. Spitler and Duncan A. Forbes 



Mieske S., Jordan A., Cote P., Kissler-Patig M., Peng 

E. W., Ferrarese L., Blakeslee J. P., Mei S., Merritt D., 

Tonry J. L., West M. J., 2006, ApJ, 653, 193 
Montuori M., Di Matteo P., Lehnert M. D., Combes P., 

Semelin B., 2010, A&A, 518, A56+ 
Ouchi M., Shimasaku K., Okamura S., Furusawa H., 

Kashikawa N., Ota K., Doi M., Hamabe M., Kimura 

M., Komiyama Y., Miyazaki M., Miyazaki S., Nakata F., 

Sekiguchi M., Yagi M., Yasuda N., 2004, ApJ, 611, 660 
Owers M. S., Randall S. W., Nulsen P. E. J., Couch W. J., 

David L. P., Kempner J. C, 2011, ApJ, 728, 27 
Peng E. W., Jordan A., Cote P., Blakeslee J. P., Ferrarese 

L., Mei S., West M. J., Merritt D., Milosavljevic M., Tonry 

J. L., 2006, ApJ, 639, 95 
Peng E. W., Jordan A., Cote P., Takamiya M., West M. J., 

Blakeslee J. P., Chen C, Ferrarese L., Mei S., Tonry J. L., 

West A. A., 2008, ApJ, 681, 197 
Puzia T. H., Zepf S. E., Kissler-Patig M., Hilker M., Minniti 

D., Goudfrooij P., 2002, A&A, 391, 453 
Romanowsky A. J., Strader J., Spitler L. R., Johnson R., 

Brodie J. P., Forbes D. A., Ponman T., 2009, AJ, 137, 

4956 

Schlegel D. J., Finkbeiner D. P., Davis M., 1998, ApJ, 500, 
525 

Schweizer F., Seitzer P., 1998, AJ, 116, 2206 

Silverman B. W., 1986, Density estimation for statistics 
and data analysis 

Spitler L. R., 2010, MNRAS, 406, 1125 

Spitler L. R., Forbes D. A., Strader J., Brodie J. P., Gal- 
lagher J. S., 2008, MNRAS, 385, 361 

Spitler L. R., Larsen S. S., Strader J., Brodie J. P., Forbes 
D. A., Beasley M. A., 2006, AJ, 132, 1593 

Strader J., Beasley M. A., Brodie J. P., 2007, AJ, 133, 2015 

Strader J., Brodie J. P., Spitler L., Beasley M. A., 2006, 
AJ, 132, 2333 

Tal T., van Dokkum P. G., Nelan J., Bezanson R., 2009, 
AJ, 138, 1417 

van den Bosch R. C. E., van de Ven G., Verolme E. K., 
Cappellari M., de Zeeuw P. T., 2008, MNRAS, 385, 647 

Woodley K. A., Harris W. E., Puzia T. H., Gomez M., 
Harris G. L. H., Geisler D., 2010, ApJ, 708, 1335 

Yagi M., Kashikawa N., Sekiguchi M., Doi M., Yasuda N., 
Shimasaku K., Okamura S., 2002, AJ, 123, 66 

Yoon S., Yi S. K., Lee Y., 2006, Science, 311, 1129 

Zepf S. E., Ashman K. M., 1993, MNRAS, 264, 611 

Zinn R., 1985, ApJ, 293, 424 



APPENDIX A: STATISTICAL TESTS FOR 
TRIMODALITY 

Al Kaye's Mixture Model algorithm 

The most common tool to determine the statistical sig- 
nificance of GC bimodality is the Kaye's Mixture Model 
(KMM) code ( [Ashman et al.|1994 1, designed to calculate the 
'p' statistic for the significance of a homoscedastic bimodal 
Gaussian distribution being preferred over a unimodal Gaus- 
sian distribution. KMM is also capable of determining the 



VCC number 
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M87 


NGC 4365 
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Table A2. KS test results comparing galaxies in the ACS Virgo 
Cluster Survey. The 'p' statistic shown is an inverse measure of 
the probability that GCs in the different galaxies are drawn from 
different underlying colour distributions. Only M49 and M87 are 
inconsistent with being drawn from different colour distributions 
and hence could be said to have the same underlying colour dis- 
tribution. Other galaxy comparisons show a probability greater 
than 99 per cent of being drawn from different underlying colour 
distributions. 



significance of heteroscedastic multimodal distributions be- 
ing preferred over a unimodal Gaussian distribution but the 
'p' statistic in these cases is an approximation. Applied to 
our GC colour distribution for NGC 4365 both bimodal and 
trimodal distributions result in a 'p' statistic less than 10"'' 
(a probability larger than 99.99 per cent that the multi- 
modal distributions are preferred over the unimodal distri- 
bution). The traditional KMM code is however not capable 
of determining the probability that a trimodal distribution 
is preferred over a bimodal distribution. 

Using an extension to the traditional KMM code, writ- 
ten by Dr Matthew Owers ( Owers et aI.|20lT |, we were able 
to assess a probability for three Gaussian modes over two 
Gaussian modes bearing in mind that the code was untested 
for application to GC colour distributions. We attempted to 
test the code at the same time as using it for our analysis 
by including checks on how the code performed with varying 
GC input numbers, analysis of other galaxies' GC distribu- 
tions and redundant calculations for the probability of four 
modes over three modes and bimodality over unimodality. 
The code showed variation of results with input GC number 
and gave results for other galaxies inconsistent with results 
in literature. Comparing results for NGC 4365 with NGC 



1407 (data from Spitler et al. 2011 in prep, also used in Ro- 
imanowsky et al.|2009| [Foster et al.|2010| [Fbrbes et al.|2011 \ 



the extended KMM test finds that both are trimodal (NGC 
4365 with a probability of 99.7 per cent and NGC 1407 with 
a probability larger than 99.9 per cent) but there is little ev- 
idence in the literature for NGC 1407 being trimodal. When 



the number of GCs is halved for NGC 4365 (data from Peng 



et al. 2006 1 the probability of trimodality decreases below 
95 per cent (to 91.8 per cent), even though the smaller sam- 
ple is more central and we expect trimodality to be more 
prominent there. 



A2 The Kolmogorov-Smirnov test 

The Kolmogorov-Smirnov (KS) test shows clearly that most 
giant elliptical galaxies have GC colour distributions that 
are significantly different from each other (subpopulation 
peak colours and widths vary between galaxies, see Table 
A2 I (using GC catalogues from Peng et al.|2006 1. A sample 
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Comparing galaxies with Subaru/S-Cam photometry and g' — i' colours where i' ^ 23 
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Table AI. Results from the extended KMM test for testing the significance with which m + 1 modes are preferred over m modes in 
colour. For each mode the mean value (fi), width (cr) and number of objects assigned (n) are shown and between each modality the 'p' 
statistic is shown. 
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Table A3. KS test results comparing simulated bimodal and tri- 
modal distributions with NGC 4365's colour distribution. The 'p' 
statistic shown is an inverse measure of the probability that the 
simulated distributions are drawn from different colour distribu- 
tions to the actual ones. 



ble |A3| shows that for all comparisons the KS test shows 
probabilities of almost 100 per cent that the simulated and 
actual distributions are significantly different. Neither sim- 
ulated distribution is a good enough description of the ob- 
served distribution, likely because they cannot take into ac- 
count the skewness and peakiness (kurtosis) of the observed 
distributions. The colour modes in GC systems are not per- 
fectly Gaussian and therefore this method of analysis cannot 
distinguish between bimodal and trimodal Gaussian descrip- 
tions of the colour distribution to determine which is more 
likely. 



of galaxies with reasonably large GC numbers were chosen to 
compare to each other and most galaxies show a probability 
of larger than 99 per cent of being drawn from a colour 
distribution different to any other galaxy in the sample. Two 
very large galaxies (M87 and M49) are consistent with being 
drawn from the same underlying colour distribution. It is 
clear that each galaxy's GC system is formed independently 
of any other GC system but that similar galaxy properties 
can result in similar GC systems. 

We also used this test to compare colour distributions 
simulated from bimodal and trimodal descriptions of the 
data (as output by KMM) with the actual GC colour dis- 
tribution and found that the observed GC distribution is 
significantly different from both simulated distributions. Ta- 



A3 Chi-Squared minimization 

Comparing the bimodal and trimodal results for a reduced 
Chi-Squared (x^) minimization is dependent on binning the 
data in colour. We have chosen to do this binning via an 
Epanechnikov kernel smoothing of the data rather than an 
histogram. With different colour bin sizes the statistical re- 
sults will also be different. We have chosen bin sizes for the 
Epanechnikov kernel that are equal to the colour errors in 
the data. 

Using the entire radial range of our data we can achieve 
the best number statistics but the results are not conclusive. 
For the ACS GCs, where we have better colour resolution, 
the best-fitting trimodal Gaussian distribution has a slightly 
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Figure Al. Results of bimodal and trimodal fits to the Epanech- 
nikov kernel smoothing of the HST/ACS GC colour distribution. 
The distribution is plotted in grey, showing poissonian errors and 
the individual Gaussians as well as the summation of all three 

arc ovcrplottcd for both the bimodal (dashed lines) and trimodal 
(dotted lines) fits. The reduced value is 1.40 for the bimodal 
case and 1.26 for the trimodal case. 



APPENDIX B: SURFACE DENSITY DATA 




Figure A2. Results of bimodal and trimodal fits to the Epanech- 
nikov kernel smoothing of the Subaru/S-Cam GC colour distri- 
bution. The distribution is plotted in grey, showing poissonian 
errors and the individual Gaussians as well as the summation of 
all three are overplotted for both the bimodal (dashed lines) and 
trimodal (dotted lines) fits. The reduced value is 2.57 for the 
bimodal case and 1.07 for the trimodal case. 



smaller reduced value than the best-fitting bimodal dis- 
tribution but the best fit third mode is very close to blue 
modes and contains many fewer objects. This 'third mode' 
is not best described as a third GC subpopulation but as an 
indication that the two modes of the underlying GC distri- 
bution are not pure Gaussian but probably Gaussians with 
both skewness and kurtosis. We do see a significantly smaller 
reduced value and a central third mode with a large pro- 
portion of the GCs in the S-Cam catalogue. 



© 0000 RAS, MNRAS 000, 000-000 



Wide-field imaging of NGC 4-365's globular cluster system: The third subpopulation revisited 27 



GC candidate radial surface density obtained from HST/ACS photometry 
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GC candidate radial surface density obtained from Subaru/S-Cam photometry 
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Table Bl. GC radial surface density data for the total GC system and eaich subpopulation. We show surface density with different 
radial binning (columns typeset in bold) for the total GC system, the blue and red subpopulations and the green subpopulation. 
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